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ABSTRACT 

An outburst of more than 80 individual bursts, similar to those seen from Soft Gamma Repeaters (SGRs), was 
detected from the anomalous X-ray pulsar (AXP) IE 2259+586 in 2002 June. Coincident with this burst activity 
were gross changes in the pulsed flux, persistent flux, energy spectrum, pulse profile, and spin-down of the 
underlying X-ray source. We present Rossi X-Ray Timing Explorer and X-Ray Multi-Mirror Mission observations 
of IE 2259+586 that show the evolution of the aforementioned source parameters during and following this 
episode and identify recovery timescales for each. Specifically, we observe an X-ray flux increase (pulsed and 
phase-averaged) by more than an order of magnitude having two distinct components. The first component is 
linked to the burst activity and decays within ~2 days, during which the energy spectrum is considerably harder 
than during the quiescent state of the source. The second component decays over the year following the glitch 
according to a power law in time with an exponent -0.22 ± 0.01. The pulsed fraction decreased initially to 
~15% rms but recovered rapidly to the preoutburst level of ~23% within the first 3 days. The pulse profile 
changed significantly during the outburst and recovered almost fully within 2 months of the outburst A glitch of 
size Ai/max/^ = (4.24 ± 0.1 l)x 10“ 6 was observed in IE 2259+586, which preceded the observed burst activity. 
The glitch could not be well fitted with a simple partial exponential recovery. An exponential rise of ~20% of the 
frequency jump with a timescale of ~14 days results in a significantly better fit to the data; however, contam- 
ination from a systematic drift in the phase of the pulse profile cannot be excluded. A fraction of the glitch 
(~19%) was recovered in a quasi-exponential manner having a recovery timescale of ~16 days. The long-term 
postglitch spin-down rate decreased in magnitude relative to the preglitch value. The changes in the source 
properties of IE 2259+586 during its 2002 outburst are shown to be qualitatively similar to changes seen during 
or following burst activity in two SGRs, thus further solidifying the common nature of SGRs and AXPs as 
magnetars. The changes in persistent emission properties of IE 2259+586 suggest that the star underwent a 
plastic deformation of the crust that simultaneously impacted the superfluid interior (crustal and possibly core 
superfluid) and the magnetosphere. Finally, the changes in persistent emission properties coincident with burst 
activity in IE 2259+586 enabled us to infer previous burst-active episodes from this and other AXPs. The 
nondetection of these outbursts by all-sky gamma-ray instruments suggests that the number of active magnetar 
candidates in our Galaxy is larger than previously thought 

Subject headings : pulsars: general — stars: individual (IE 2259+586) — stars: neutron — X-rays: bursts 


1. INTRODUCTION 

Anomalous X-ray pulsars (AXPs) and soft gamma repeaters 
(SGRs) are two intriguing classes of isolated neutron stars, 
very likely magnetars, whose bright X-ray emission is pow- 
ered by the decay of their strong magnetic fields. When the 
common nature of AXPs and SGRs was first proposed by 
Thompson & Duncan (1996) with both being magnetars, the 
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observational evidence linking them was tenuous. A major 
advance in connecting these two classes came when slow 
pulsations and rapid spin-down, defining characteristics of 
AXPs, were discovered from SGR 1806—20 (Kouveliotou 
et al. 1998). Since then, further observational similarities 
among AXPs and SGRs have been established. For example, 
Marsden & White (2001) performed a systematic analysis of 
SGR and AXP spectral data that showed the two classes form 
a continuum in spectral hardness versus spin-down rate in 
which the SGRs have harder spectra and faster spin-down 
rates than the spectrally softer, slower braking AXPs. Simi- 
larly, the timing noise strength in SGRs and AXPs appears 
to be correlated with spin-down rate (Woods et al. 2000; 
Gavriil & Kaspi 2002). For a review of SGRs and AXPs, see 
Woods (2004) and Mereghetti et al. (2002), respectively. 
Until recently, the most prominent characteristic distin- 
guishing the two groups was the emission of soft bright (up to 
10 44 ergs s” 1 ) bursts of soft gamma-rays from the SGRs and 
not the AXPs. In fact, it was this extraordinary property of 
SGRs that led to doubt within the community that AXPs and 
SGRs were of the same nature. This uncertainty was removed 
when SGR-like bursts were recently discovered from at least 
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one AXP (Kaspi et al. 2003) and probably one other (Gavriil, 
Kaspi, & Woods 2002). 

SGR-like bursts from the direction of an AXP were first 
detected from the source IE 1048. 1 -5937 (Gavriil et al. 2002). 
A single weak SGR-like burst was detected during each of two 
Rossi X-Ray Timing Explorer (RXTE) Proportional Counter 
Array (PCA) monitoring observations of this AXP separated 
by 2 weeks. However, the identification of IE 1048.1-5937 
as the burst source could not be made unambiguously because 
of the lack of imaging capability with the PCA. Interestingly, 
the quiescent properties of this AXP (e.g., energy spectrum, 
pulse profile, and timing noise) most closely resemble those of 
the SGRs (Kaspi et al. 2001), making this AXP the most SGR- 
like of its class. 

The second detection of SGR-like bursts was recorded from 
the AXP IE 2259+586 on 2002 June 18. This is the least 
SGR-like of the AXPs in terms of its quiescent source prop- 
erties. Unlike the two weak bursts observed earlier from IE 
1048.1-5937, this AXP showed a major SGR-like outburst or 
collection of bursts (Kaspi et al. 2003). In total, more than 80 
bursts were detected within 3 hr of observing time. A detailed 
study of these bursts will be presented in a companion paper 
(Gavriil, Kaspi, & Woods 2003). In addition to emitting 
these hard X-ray bursts, several parameters of the persistent 
source changed in conjunction with this outburst, thereby 
confirming the AXP as the source of the burst emission (Kaspi 
et al. 2003). 

Here we present X-Ray Multi-Mirror Mission (XMM- 
Newton) and RXTE PCA observations of the persistent X-ray 
flux from IE 2259+586 before, during, and after the 2002 June 
outburst. We quantify the changes of the spectral and temporal 
properties of the X-ray source, as well as the timescales for 
their recovery. We compare the changes observed in this AXP 
with dynamic behavior seen in the persistent emission of SGR 
1900+14 (Woods et al. 2001) and SGR 1627-41 (Kouveliotou 
et al. 2003) during episodes of intense burst activity. Finally, 
we present a possible explanation for the observed behavior 
in IE 2259+586 within the context of the magnetar model. 

2. XMM-NEWTON OBSERVATIONS 

The results presented in this section were obtained from 
observations of IE 2259+586 with the telescopes aboard the 
XMM-Newton observatory (Jansen et al. 2001). This observa- 
tory is composed of three co-aligned X-ray telescopes. The 
focal plane instruments are one EPIC pn camera (Striider et al. 
2001) and two EPIC MOS cameras (Turner et al. 2001). All 
instruments are sensitive to X-rays between 0.2 and 15.0 keV. 
The pn camera has an effective area of ~1400 cm 2 at 1.5 keV, 
while the MOS cameras each have areas of ~500 cm 2 at 
1 .5 keV. The focused X-ray beam for the telescopes serving 


each MOS camera is split between the focal plane instrument 
and the Reflection Grating Spectrometer (RGS; den Herder et al. 
2001 ). 

The XMM-Newton observatory observed IE 2259+586 5 
times during 2002. Three pointings were focused on different 
portions of CTB 109, the supernova remnant (SNR) sur- 
rounding IE 2259+586. The primary scientific objective of the 
remaining two observations was the central point source. In 
each of the three CTB 109 pointings, the AXP is within the 
field of view at large off-axis angles (10'- 13')- Results from 
all five observations on the central source will be presented 
here. Analysis of the data from the remnant will be presented 
elsewhere (M. Sasaki et al. 2004, in preparation). 

The first XMM-Newton observation of IE 2259+586 was 
one of the remnant pointings, carried out on 2002 January 22. 
The second was centered on the AXP and was performed on 
2002 June 11, fortuitously 1 week prior to the outburst of IE 
2259+586. Following the outburst, a target of opportunity (ToO) 
was declared and the source was reobserved 3 days later on 2002 
June 21. The remaining two observations of CTB 109 were 
carried out on 2002 July 9. See Table 1 for further details on these 
observations. Hereafter, each XMM-Newton observation is re- 
ferred to by the identifier label assigned to it in Table 1. 

The two MOS cameras were operated in full-firame mode 
for all but one observation (Obs3, Small Window for MOS1) 
to study the SNR. The frame time for the full-frame data mode 
is 2.6 s, which causes severe pile-up for the AXP (£50%). 
Therefore, these data are not considered further here. 

The pn camera was operated in small-window mode for 
both Obs2 and Obs3, the two observations closely bracketing 
the 2002 June outburst in time. The frame time for the pn 
camera in this mode is ~5.86 ms, allowing detailed study of 
the pulsed emission and a search for burst emission. For the 
on-axis count rate of the AXP, the dead-time fraction is 30%, 
but the pile-up fraction is only 0.03%. The three off-axis 
pointings (Obsl, Obs4, and Obs5) were carried out in ex- 
tended full-frame mode, which has a time resolution of 
200 ms. The high off-axis angles for the AXP reduced the count 
rate by a factor of 2-3. The different data mode and reduction 
in count rate lessened the dead time to 2%, but the pile-up 
fraction increased to 3%. 

The RGS data from Obs2 and Obs3 were acquired in spec- 
troscopy mode with 5.7 s time resolution and excellent energy 
resolution (A E/E = 10-20). The coarse time resolution al- 
lowed for only phase-averaged spectral analysis (see § 2.3). 

All data were processed using the XMM-Newton Science 
Analysis System 8 ^ (S AS) version 5.4.1. For the pn data, the 


8 See http://xmm. vilspa.esa.es/external/xmm_sw_cal/sas.shtml. 


TABLE 1 

XMM-Newton Observation Log for IE 2259+586 


Name 

XMM-Newton ObsED 

Relation to Burst Activity 

Date* 

(MJDTDB) 

pn Exposure 11 
(ks) 

pn Data Mode c 

AXP Offset* 1 
(arcmin) 

Obsl 

0057540101 

Before 

52,296.791 

10.7 

Extended 

11.2 

Obs2 

0038140101 

Before 

52,436.546 

24.9 

Small 

1.7 

Obs3 

0155350301 

After 

52,446.449 

18.5 

Small 

1.7 

Obs4 

0057540201 

After 

52,464.369 

12.3 

Extended 

13.0 

Obs5 

0057540301 

After 

52,464.606 

12.4 

Extended 

9.7 


• Start time of each observation. Note that the outburst began on date 52,443.65 MJD. 

b Exposure times quoted reflect on-source times after filtering of flares, etc. for spectral and temporal analysis. 
c Extended full-frame mode or small-window mode. 

d The angular distance of IE 2259+586 from the optical axis of the EPIC mirrors. 
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TABLE 2 

Pulse Frequencies and Pulsed Fractions op IE 2259+586 Measured Using XMM-Newton pn Data 


Observation 

Epoch 

(MJDTDB) 

Frequency* 

(Hz) 



Pulsed Fraction 



0.3- 1.0 keV 

1 .0-2.0 keV 

2.0-5.0 keV 

5.0-12.0 keV 

2.0-10.0 keV 

Obsl 

52,296.79 1 -52,296.932 

0.14328688 (55) 

0.169 (15) 

0.195 (6) 

0.230 (10) 

0.314 (73) 

0.232 (10) 

Obs2 

52,436.546-52,436.983 

0.14328705 (13) 

0.215 (6) 

0.225 (3) 

0.234 (4) 

0.284 (23) 

0.234 (4) 

Obs3 

52,446.449-52,446.755 

0.14328746 (9) 

0.168 (5) 

0.200 (2) 

0.223 (3) 

0.189 (13) 

0.220 (3) 

Obs4 b 

52,464.369-52,464.534 

0.14328771 (13) 

0.166 (15) 

0.180 (7) 

0.225 (11) 

0.294 (50) 

0.230 (10) 

Obs5 b 

52,464.606-52,464.761 

0.14328771 (13) 

0.139 (13) 

0.168 (6) 

0.225 (8) 

0.339 (38) 

0.230 (8) 


* Numbers in parenthses indicate the 1 a uncertainty in the least significant digits of the frequency. 
b Frequency derived from combined data set of Obs4 and Obs5 due to their close proximity to one another in time. 


script EPCHAIN was run on the observation data files. This 
script processes the data for use in higher end analysis tools. A 
light curve of the full field of view (FOV) excluding the bright 
central source was constructed (0.5-10.0 keV) and inspected 
for times of high background. We chose a threshold of 2 times 
the nominal background to define regions of high background. 
Periods of high background constituted 0%-33% of each data 
set and were eliminated from further analysis. 

2.1. Burst Search 

We have used the data from the pn camera to search for 
burst emission from IE 2259+586 during the XMM-Newton 
observations. By using the SAS tool E VS ELECT, events from 
within 10" and 12?5 radius circles around the position of the 
AXP were extracted for the on-axis and off-axis (remnant) 
pointings, respectively. Following the standard filtering pro- 
cedures for XMM-Newton pn data, grades >12 and flag values 
equal to 0 were retained. Next, we filtered the event list for 
energy between 0.5 and 12.0 keV to optimize the signal-to- 
noise ratio for burst emission. 

Light curves were constructed for each of the XMM-Newton 
observations at the nominal time resolution (6 ms for Obs2 
and Obs3 and 200 ms for the others), 0. 1 s resolution (Obs2 and 
Obs3), 1 s, and 10 s resolution. No significant burst emission 
was detected on these timescales during any of these obser- 
vations. This is consistent with the absence of bursts in the 
more densely sampled RXTE data (Gavriil et al. 2003) that 
bracket these XMM-Newton observations. Bursts were re- 
corded from IE 2259+586 only on 2002 June 18. 

2.2. Pulse-Timing Analysis 

Using the same procedures as those described above for the 
burst search, source event lists from the pn data were extracted 
to study the pulsed emission from IE 2259+586. To optimize 
the signal-to-noise ratio for the pulsed emission, the event list 
was filtered for energy between 0.5 and 7.0 keV. 

There is a known problem with the time tags of XMM- 
Newton pn data in that there can be sudden jumps in the 
photon arrival times of integral second values (W. Becker 
2003, private communication). It is still not understood at 
what stage of the processing these “time jumps” distort the 
time tags, only that they do occur. We searched for and iden- 
tified four time jumps within two of the five data sets (Obs 2 
and 3). All time jumps were identified in the two pointings 
at which the pn camera was operated in small-window 
mode. To find these time jumps we calculated the modulus of 
each event time stamp with the mean frame time of the cor- 
responding data set. The frame time varied with data mode 


(small window vs. extended full frame), as well as between and 
within observations having identical data modes. The latter 
effect is believed to be due to temperature variations in the on- 
board electronics {XMM-Newton Helpdesk 2003, private 
communication). 

The time jumps in each data set were obvious after plotting 
the frame time residuals versus the time of each event Time 
jumps manifested themselves as discontinuities in the frame 
time residuals plot These discontinuities were identified and 
corrected for by adding or subtracting an integral number of 
seconds to the data following the jump until the frame time 
residuals matched precisely across the boundary. The time 
corrections we applied to the data were verified by the pulsar 
data. In the uncorrected pulsar data set, phase jumps that were 
consistent with being equal to the time corrections required by 
the frame time inconsistencies were detected. By using the 
time-corrected data set, the phase jumps in the pulsar data 
disappeared. Therefore, we are secure about the relative tim- 
ing of the corrected XMM-Newton pn data set. However, the 
absolute timing of these events will require a better under- 
standing of the origin of the time jump problem. 

The corrected time tags were next converted to the solar 
system barycenter by using the SAS tool BARYCORR. As- 
suming that the remaining uncertainty in the corrected time 
tags is a small integral number of seconds, the propagated 
error in the barycenter correction applied to the data is in- 
significant with respect to the precision with which we can 
measure the pulse frequency. 

For the two point-source observations in small-window 
mode (Obs2 and Obs3), the data were binned at twice the 
nominal frame time of each observation, and a fast Fourier 
transform was constructed of the light curves (0.5-7.0 keV). 
The only significant power detected between 4 x 10 -4 Hz 
and the Nyquist frequency (44 Hz) is from the first seven or 
eight harmonics of the pulsar. The 3 a upper limit to the 
mean fractional power in the frequency range 0.001-1.0 Hz 
(minus the pulsar frequency and harmonics) is 3 x 10 -3 and 
2 x 10~ 3 (rms/mean) 2 Hz -1 , for Obs2 and Obs3, respectively. 
These limits are orders of magnitude lower than typical 
broadband noise power levels seen in accreting X-ray pul- 
sars, even so-called quiet accreting pulsars such as 4U 1626— 
67 (Owens, Oosterbroek, & Parmar 1997; Shinoda et al. 
1990). 

The precise pulse frequencies for each pointing were de- 
termined using a phase-folding technique. For each data set, a 
template pulse profile was constructed by folding the data on 
the frequency as determined from the peak power of the 
fundamental frequency in the power density spectrum. Next, 
the data set was split into wlO 3 s segments, and each segment 
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was folded with the same frequency. The pulse profiles from 
each segment were cross-correlated with the template pulse 
profile, and a phase offset was measured. The resulting phase 
differences for each segment were fitted to a line, and the slope 
was used to correct the frequency. By using the refined fre- 
quency, a new template pulse profile was constructed, and the 
procedure for refining the pulse frequency was repeated. The 
epochs and frequencies measured for each XMM-Newton ob- 
servation are listed in Table 2. Since two of the CTB 109 
observations were performed sequentially, we determined a 
single frequency for the combined observation. The spin fre- 
quencies measured in the XMM-Newton observations are 
consistent with the much more precise spin ephemeris mea- 
sured using the RXTE data (see § 3.1). 

Using the measured frequencies, the background-subtracted 
pn data were folded over different energy bands. The back- 
ground was estimated by measuring the count rate for the 
different energy bands in a circular extraction region close to 
the pulsar and within CTB 109 but away from bright knots 
within the remnant. The normalized pulse profiles for different 
energy bands are shown in Figure 1. Clearly, the pulse profile 
changed significantly directly after the outburst. Subsequent 
X-ray observations revealed a gradual recovery to the pre- 
outburst pulse shape (see § 3.4). 

The rms pulsed fractions were calculated for each obser- 
vation over several energy bands. Following van der Klis 
(1989), the pulsed fractions were calculated using the first 
seven harmonics of the Fourier representation of the pulse 


profile according to the equation 


F rms — 


y ■><** + #* ( <J a t +a %) (j) 




where 


1 n x n 

at= jvE r, ' cos 2?r< M- t 3 k= jj'Yi n sin 27r ^+ 

V i-\ V i=l 

1 N l N 

< c< * 2 2 nfrk, <4 = jyl XX sin 2 2n<f>ik. 


Here Frms is the pulsed fraction, k refers to the harmonic 
number, i refers to the phase bin, N is the total number of 
phase bins, 0, is the phase, r, is the count rate in the z'th phase 
bin, and o Xi is the uncertainty in the count rate of the zth phase 
bin. Note that equation (1) is insensitive to reducing the num- 
ber of harmonics used in calculating F ^ for low signal-to- 
noise ratio pulse profiles since the statistical noise is sub- 
tracted from the total variance. Seven harmonics were chosen 
for N to encompass all statistically significant source power in 
the highest signal-to-noise ratio pulse profiles. The pulsed 
fractions and uncertainties are given in Table 2. 

Directly following the outburst (Obs3), the pulsed fraction 
at all energies dropped relative to the preoutburst level (Obs2). 


Obsl 0bs2 Obs3 0bs4 0bs5 



0 0.5 1 1.5 2 0 0.5 1 1.5 2 0 0.5 1 1.5 2 0 0.5 1 1.5 2 0 0.5 1 1.5 2 

Phase (cycles) 


Fig. 1. — Folded-pulse profile of IE 2259+586 as observed with the XMM-Newton pn camera as a function of energy for the all five observations. Time increases 
from left to right. The burst activity is recorded between Obs2 and Obs3. See Table 1 for exact times of observations. 
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The largest change was seen in the hard band (5.0-12.0 keV). 
During the XMM-Newton observations 3 weeks after die 
outburst (Obs4 and Obs5), the broadband pulsed fraction (2- 
10 keV) had recovered to its preoutburst value, although the 
pulse profile was still significantly different The time evolu- 
tion of the pulsed fraction is presented in § 3.3, in which 
additional pulsed fraction measurements made with the RXTE 
PCA are reported. 

In the observation directly before the outburst (Obs2), the 
pulsed fraction does not vary strongly with photon energy. 
The energy dependence of the pulsed fraction is most prom- 
inent in the observations following the outburst for which the 
pulsed fraction increased significantly with energy. Interest- 
ingly, the pulsed fraction in the soft band (0.3-1.0 keV) is 
higher 1 week before the outburst than at any of the other 
XMM-Newton epochs before or after. 

2.3. Spectroscopy 

The pn data were used for spectral analysis of IE 2259+586 
because of their excellent signal-to-noise ratio and negligible 
pile-up. For the two observations centered on the point source, 
events from within a 10" radius circle around the position of 
the AXP were extracted using the SAS tool EVSELECT. 
Following the standard procedures for XMM-Newton pn data, 
these events were filtered for grades >4, and flag values 
equal to 0 were retained for the source spectrum. Similarly, a 
background spectrum was extracted using the same filtering 
criteria from nearby circular regions of 10" radius for each ob- 
servation identical to the regions used for the pulsed-fraction 
analysis. The total number of source counts accumulated for 
each spectrum was 259,000 and 373,000 for the preburst 
(Obs2) and postburst (Obs3) observations, respectively. 

The three observations of CTB 109 included the AXP but 
at large off-axis angles (10 / -13')- At off-axis positions of 
~10', the point-spread function of the EPIC mirrors is sig- 
nificantly broader than on-axis. Because of this effect, a larger 
source extraction radius (12*5) was used. Background spectra 
were extracted from circular regions of radius ~22" with 
centers at off-axis positions similar to the source extraction 
region. The larger extraction radii were chosen to increase the 
number of background counts and improve the accuracy of 
the background subtraction. Furthermore, the background 
spectra were extracted from the same chip. Because of the off- 
axis positions of the AXP and the lower exposures of the CTB 
1 09 observations, the total number of source counts extracted 
in each spectrum was much lower, between 38,000 and 
56,000. 


Each spectrum was grouped such that there was a minimum 
of 25 counts per bin. The spectra were fitted individually using 
XSPEC 9 version 1 1 .2.0. Because of the high column along the 
line of sight to this source, the observed counts in pulse height 
amplitude (PHA) channels corresponding to E < 0.6 keV are 
dominated by the low PHA tails of events whose true energies 
are above 1 keV. We, therefore, restrict our spectral fits to the 
data in the range 0.6-12.0 keV. Fits to single-component 
models (blackbody, power law, and bremsstrahlung) modified 
by interstellar absorption did not yield statistically acceptable 
fits. The resulting reduced \ 2 values to the single-component 
models were £4. We next tried the standard AXP two- 
component spectral model of a blackbody plus a power law 
modified by interstellar absorption. We obtained good fits to 
all five data sets. The results of these fits are given in Table 3. 
The X-ray spectrum of IE 2259+586 during Obs2 and the 
residuals from the best fit are shown in Figure 2. Formally, the 
fit to the Obs3 spectrum is not statistically acceptable. How- 
ever, the residuals between 0.8 and 3.0 keV constitute the 
majority of the total \ 2 and are at the few percent level. These 
residuals are very likely due to uncertainties in the instru- 
mental response (Briel et al. 2004). Introducing a 3% sys- 
tematic error in the spectral model reduces xt to unity. The 
errors quoted in Table 3 are inflated according the systematic 
error. No narrow-line features (absorption or emission) are 
evident in any of the phase-averaged spectra. Within the en- 
ergy range 0.9-2.0 keV, the 90% confidence limit on a narrow- 
line feature on the order of the pn response function (~50 eV 
FWHM) is ~10 eV for the equivalent width. Between 2.0 and 
7.0 keV, the response function increases up to ~130 eV and 
the line limit increases from ~10 to ~70 eV. 

We next included the RGS data in our spectral analysis. The 
RGS spectra for Obs2 and Obs3 were extracted using standard 
processing techniques for a point source. The data were 
grouped to 25 counts per spectral bin and ported into XSPEC 
for simultaneous fitting with the pn data. Each observation 
was fitted independently, again to the blackbody plus power 
law model, and the measured spectral parameters were con- 
sistent with those obtained using the pn data alone. The su- 
perior energy resolution of the RGS data combined with the 
high signal-to-noise ratio data from the pn camera allow us to 
put even more constraining limits on the presence of narrow- 
line features. Within the energy range 0.8-1.75 keV, the 90% 
confidence limit on a narrow-line feature (~3— 10 eV) is ~7 eV 
for the equivalent width. 

9 See http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec. 


TABLE 3 

Phase-averaged Spectral Fit Parameters op IE 2259+586 from XMM-Newton pn Data 


Parameter* 

Obsl 

Obs2 

Obs3 

Obs4 

Obs5 

Nfi (10 22 cm~ 2 ) 

1.096 (24) 

1.098 (12) 

1.035 (8) 

0.953 (24) 

0.937 (20) 

kT (ke V) 

0.488 (14) 

0.411 (4) 

0.517(5) 

0.537 (19) 

0.548 (13) 

r 

4.04 (8) 

4.10 (3) 

3.59 (2) 

3.62 (7) 

3.58 (6) 

Flux (10" 11 ergs cm" 2 s" 1 ) 1 * 

1.24 (4) 

1.30 (2) 

3.47 (3) 

2.01 (5) 

2.11 (5) 

Unabsorbed flux (10" H ergs cm -2 s _1 ) c 

153 (4) 

1.63 (2) 

4.17 (4) 

2.37 (6) 

2.49 (5) 

PL/BB ratio 11 

0.70 (6) 

0.43 (2) 

0.74 (2) 

0.93 (9) 

0.79 (6) 

X 2 /dof 

619/578 

922/851 

1332/1094 

606/540 

599/638 


* Numbers in parentheses indicate the 1 o uncertainty in the least significant digits of the spectral parameter. Note that these uncertainties reflect the 1 a 
error for a reduced x 2 of unity. 

b Observed flux from both spectral components 2-10 keV. 
c Unabsorbed flux from both spectral components 2-10 keV. 

d The ratio of the power-law flux (2—10 keV) to the bolometric blackbody flux (corrected for absorption). 
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Fig. 2. — Phase-averaged X-ray spectrum of IE 2259+586 as observed with 
the XMM-Newton pn camera during the observation 1 week prior to the 2002 
June outburst (Obs2). Bottom , Ratio of the data to the best-fit model (see 
Table 3) folded through the instrumental response. 

We used the RGS data to confirm whether the fit value of 
A H was affected by the steep power-law component of the 
spectral model. Specifically, using just the pn data, it is diffi- 
cult to distinguish between a steep power law that is strongly 
absorbed and a fiat or inverted power law (e.g., blackbody) 
that is observed through modest absorption. The Ne-K edge at 
0.87 keV is readily detected in the RGS data. Freeing both the 
TVh and the Ne abundance relative to solar, a Ne , and using the 
TBVARABS model in XSPEC, we find that a^ c = 1.27 ± 
0.24 (at 90% confidence for 1 degree of freedom [dof]) and 
that the best-fit value of N^ is very close to the best-fit value 
when UNe i s fixed at 1 • Thus, unless the true Ne abundance is 
significantly different from the solar value, the fitted A H seems 
accurate (see Table 3). 

Clearly, the energy spectrum of IE 2259+586 hardened 
following the 2002 June outburst. Between the XMM-Newton 
observation 1 week prior to the outburst and 3 days following, 
the photon index became significantly flatter and the black- 
body temperature rose. However, the blackbody temperature 
during Obsl, several months before the outburst, was signif- 
icantly higher than the Obs2 temperature and only marginally 
lower than the postoutburst temperature, albeit with a con- 
siderably lower flux. This might suggest that the low tem- 
perature measured during Obs2 was an indicator of the 
impending outburst. However, the blackbody temperature 
measured using the Chandra X-Ray Observatory in 2000 
January was 0.412 db 0.006 keV (Patel et al. 2001), and there 
was no glitch detected from IE 2259+586 at this time (Gavriil 
& Kaspi 2002). Moreover, Patel et al. (2001) measure a photon 
index of 3.6 ±0.1, consistent with the postoutburst value. 
Other previous observations of this AXP have shown similar 
variance in these spectral parameters (e.g., Parmar et al. 1998; 
Marsden & White 2001), in addition to the power law- 
to-blackbody (PL/BB) flux ratio (Marsden & White 2001). 
Thus, it appears that IE 2259+586 undergoes significant 
spectral changes independent of large outbursts, and the spec- 
tral parameters measured during all the XMM-Newton obser- 
vations (even those following the outburst) are within the 
historical range of these parameters. 

Using the frequencies given in the last section, we extracted 
spectra for 16 phase intervals per cycle per observation. For 
the observation 1 week prior to the outburst (Obs2), between 
10,000 and 22,000 counts were accumulated per phase interval. 


The spectrum having the least number of total counts (i.e., 
pulse minimum) was grouped such that at least 25 counts were 
contained in each data bin. The grouping for each of the 
remaining data sets was forced to be identical to the grouping 
used for pulse minimum. As with the phase-averaged spec- 
trum, we fitted a blackbody plus a power law modified by 
interstellar absorption. The 16 intervals were fitted simulta- 
neously with the column density forced to be the same for all 
intervals. The measured column density is consistent with the 
phase-averaged value. All other parameters were allowed to 
vary in the fit. The folded pulse profile, photon index, 
blackbody temperature, and ratio of power-law to blackbody 
flux are shown in Figure 3. The average values of the mea- 
sured spectral parameters are denoted by the horizontal dashed 
line shown within each panel. The flux ratio uncertainties 
are determined by propagating the flux errors in each spectral 
component. The probability that the variance in the spectral 
parameters with phase is purely statistical is 1.4 x 10”, 0.10, 
and less than 10 -16 for the photon index, blackbody temper- 
ature, and the flux ratio of the two components, respectively. 
The superb statistics of the pn data enabled us to identify both 
the photon index and flux ratio as varying significantly as a 
function of pulse phase, whereas the blackbody temperature is 
consistent with remaining constant. 



Phase (cycles) 


Fig. 3. — Phase-resolved X-ray spectrum of IE 2259+586 as observed with 
the XMM-Newton pn camera 7 days prior to the observed burst activity 
(Obs2). Top to bottom: Folded pulse profile over the energy range 0.5-7. 0 keV, 
die photon index, the blackbody temperature (kT\ and the ratio of power-law 
to blackbody flux. The horizontal dashed lines in the bottom three panels 
denote the average value for each spectral parameter. For the ratio of power- 
law to blackbody flux, the power-law flux was summed over the energy range 
2-10 keV, and the blackbody was summed over all photon energies, each 
corrected for the interstellar absorption (see text for value). 
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Fig. 4.— Same as Fig. 3, but 3 days after the observed burst activity (Obs3) 


Taken at face value, this observation favors two distinct 
components to the IE 2259+586 spectrum as has been argued 
previously (Thompson & Duncan 1996; Pema et al. 2001). 
However, the two components are likely highly correlated to 
explain the minimal variation of the pulse fraction and shape 
versus energy (Ozel, Psaltis, & Kaspi 2001). Alternatively, a 
single thermal component modified by the strong magnetic 
field (Ozel 2001) may be possible if there is substantial var- 
iation in the magnetic field across the stellar surface. 

The same analysis procedure was applied to the data directly 
following the outburst The folded profile and spectral param- 
eters versus pulse phase are shown in Figure 4. Unlike the 
observation 10 days earlier, the changes in photon index are 
only marginally significant (2.9 x 10 -3 ). Similarly, the black- 
body temperature changes are also marginally significant 
(1 x 10 -3 ). The flux ratio, on die other hand, still shows signif- 
icant variability (<10 -16 ). The phase dependence of the flux 
ratio data is markedly different from what was seen preoutburst 

The total number of counts recorded from IE 2259+586 
within the off-axis pointings (Obsl/4/5) is comparable to the 
counts recorded within only two phase bins for the on-axis 
pointings. Because of the poorer count statistics in these 
observations, phase-resolved spectroscopy was not performed 
on these data. 


3. RXTE PCA OBSERVATIONS 

The bursting behavior on 2002 June 18 from IE 2259+586 
(Fig. 5) was detected as part of an ongoing monitoring cam- 
paign of AXPs spanning the last 6 years (e.g., Kaspi, 
Chakrabarty, & Steinberger 1999; Gavriil & Kaspi 2002). 
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Fig. 5 . — RXTE PCA light curve of IE 2259+586 recorded on 2002 June 1 8. 
The light curve displays counts within 2-20 keV at 1 s time resolution. 


Follow-up ToO observations of the source with the PCA were 
immediately triggered and started as early as 1 day following 
the outburst. ToO observations continued for roughly the next 
month before regular monitoring observations resumed. Here 
we report on observations of the persistent and pulsed emis- 
sion from the AXP leading up to, during, and following the 
2002 outburst of this source. A detailed analysis of the burst 
properties is presented in a companion paper (Gavriil et al. 
2003). 

The PCA instrument aboard RXTE is a collimated (1° 
FWHM FOV) proportional counter containing a mixture of 
xenon and methane gas (Jahoda et al. 1996). It is sensitive to 
X-rays in the 2-60 keV bandpass and has a maximum effec- 
tive area of ~6500 cm 2 at 7 keV. The data from the PCA 
are read out in a number of different data acquisition modes. 
Different data modes were used depending on the analysis per- 
formed. In all the analysis, integration times including de- 
tected burst emission were eliminated from the accumulated 
spectra and/or light curves. 

With the exception of die first two ToO pointings directly 
following die outburst, GoodXenonWithPropane data were 
acquired during each observation. These data were used for 
the pulse-timing analysis described below. For the first 
two ToO observations, data modes better suited for studying 
bright bursts were employed; however, no further burst ac- 
tivity was seen. For these pointings, the event data mode 
(E_125us_64M_(M s) was used for the timing analysis. 

Starting from RXTE production level data, GoodXenon- 
WithPropane and event mode data were energy selected 
(2-10 keV) for all xenon layers and binned into light curves 
having a 0.0625 s time resolution. The time values in the light 
curve were then corrected to the solar system barycenter. For 
these processing steps, the standard prescriptions for RXTE 
PCA data analysis were followed. 10 The data were filtered 
further by eliminating times of high background and bursts 
during the 2002 June 18 observation by using custom soft- 
ware. The count rates were normalized to the number of 
proportional counter units (PCUs) on at a given time. 


10 


See http://heasarc.gsfc.nasa.gov/docs/xte/recipes/cook_book.html. 
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3.1. Pulse-Timing Analysis 

Within 2 days of the outburst, it was clear from the PCA 
data that a glitch had occurred (Kaspi et al. 2003). A new 
timing ephemeris was established and refined with continuing 
PCA observations. It gradually became apparent that accom- 
panying the sudden increase in frequency was a dramatic in- 
crease in spin-down rate (i.e., torque) by a factor of ~2. This 
torque excess decayed rapidly over the next several weeks, 
approaching the preoutburst level. Here we extend the pulse- 
timing ephemeris 15 months beyond the glitch, allowing us to 
better characterize the frequency evolution postglitch and 
quantify the recovery timescale of the torque. 

As with the XMM-Newton pn data, a phase-folding tech- 
nique was used to determine die precise ephemeris from the 
RXTE PCA data. Briefly, photon arrival times, obtained using 
data in the energy range 2-10 keV, were binned with 62.5 ms 
time resolution, reduced to the solar system barycenter, and 
then folded at the nominal pulse period. The folded pulse 
profiles for each pointing were cross-correlated in the Fourier 
domain with a template pulse profile obtained from preoutburst 
observations only, and a relative phase (i.e., time of arrival) 
was measured. Our analysis includes a total of 62 phase mea- 
surements obtained between 2000 March and 2003 September, 
with 43 measurements obtained either during or postoutburst. 

This particular analysis was complicated by the change in 
the pulse profile at the time of the outburst (see § 3.4), as the 
cross-correlation procedure assumes a fixed profile. Gross 
pulse profile changes in which the relative amplitudes of the 
two peaks were reversed resulted in a misidentification of the 
standard fiducial point by the cross-correlator for a subset of 
the phase measurements made during the outburst. This was 
accounted for by aligning the average pulse profile during the 
outburst observation with the template profile. The aligned 
outburst profile was then used as the template for phase 
measurements during the outburst. We verified this procedure 
by comparing the phase measurements made using this method 
with those calculated using the original template. We found 
that when the cross-correlator identified the appropriate peak 
in the original measurements, the revised measurements agreed 
to within the errors. We further checked that our method 
worked by verifying that the results we describe below are 
insensitive to the omission of the phase measurements made 
within 1 day of the burst activity for which the pulse profile 
changes were largest. 

We attempted to fit the spin evolution through and following 
the glitch with standard glitch models, that is, a simple jump in 
frequency, and a two-component frequency jump in which one 
part decays exponentially (eq. [2] with Av g = 0). Neither 
provided a good fit to the phase data. Figure 6 shows the fre- 
quency evolution and phase residuals following subtraction of 
die best-fit model including a sudden frequency jump and an 
additional exponential decay. The best-fit parameters imply a 
total frequency jump (Ai/) of 6.3 x 10~ 7 Hz, with a fraction 
Q = 0.23 decaying on a timescale of ~40 days. However, as 
can be seen from the figure, the residuals from this model show 
a significant systematic trend at the few percent level 
(X 2 = 447 for 54 dof). Omitting the immediate postburst data 
in which the pulse profile had substantially changed does not 
alter the result. It is possible, however, that the residuals are a 
result of low-level systematic pulse morphology variations (see 
§ 3.4). This is hard to rule out. 

However, as we show next, a model in which a substantial 
portion of the glitch is resolved in time provides a better fit to 


2000 2001 2002 2003 



Time (MJO TDB) 

Fig. 6. — Top: Frequency evolution of IE 2259+586 around the time of the 
outburst for a simple exponential recovery model. See text for details. The line 
represents the best-fit model. The circles denote frequency measurements for 
independent subsets of data. Middle: Frequency residuals of the independent 
frequency points minus the model. Bottom : Phase residuals with respect to the 
best-fit model. Clearly, there is a systematic trend in the postglitch residuals 
reflecting the inadequacy of this model. 

the data. We employed a more complex model previously 
invoked for fitting two glitches from the Crab pulsar in which 
the rise in frequency is resolved (Lyne, Pritchard, & Graham- 
Smith 1993; Wong, Backer, & Lyne 2001). We model the fre- 
quency evolution as 

v — i/ 0 (t) + Av 4- Av g {\ — 

- Ai/ rf (l - e+-‘»)/ r ') + AOt, (2) 

where i/o(t) is the frequency evolution preglitch, Av is a in- 
stantaneous frequency jump, Av g is the resolved frequency 
jump that grows exponentially on a timescale r g , Ai/j is the 
postglitch frequency drop that decays exponentially on a time- 
scale Trf, t g is the glitch epoch, and At> is the postglitch change 
in the long-term frequency derivative. 

The full ensemble of phases was fitted to the model above 
using a Levenberg-Marquardt least squares-fitting routine. 
The fit improved significantly over the simple jump in fre- 
quency and partial exponential recovery (x 2 = 131 for 53 dof). 
Because of the strong correlation between the amplitudes of the 
two exponential factors (Ai^ and AvJ) and similar associated 
timescales, our fit was mildly nonconvergent. For this reason, 
we can quote only lower limits for each amplitude. Fixing 
either one of the exponential amplitudes to a value above its 
lower limit allows the fit to converge. The key parameter is the 
difference between the two amplitudes, which determines the 
peak frequency following the glitch. The best-fit parame- 
ters for this fit (including the relationship between Av g and 
A Vd) are given in Table 4. The variance of the timescales ( r g 
and Td ) is far less than the amplitudes. We quote the for- 
mal errors for these timescales by using a fixed Av g = 
2.3 x 10~ 6 Hz. However, fixing Av g to the minimum allowed 
value yields significantly different timescales [r g = 12.8 (7) 
days and tj = 17.4 (6) days]. Going to higher values for Av g 
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TABLE 4 

Spin Parameters for IE 2259+586 from 3.2 yr of Phase-Coherent 
Timing Using RXTE PCA Data 


Parameter 

Value 

Spin frequency, i / 

0.14328703257 (21) Hz 

Spin frequency derivative, v 

-9.920 (6) x 10-' 5 Hz s" 1 

Epoch 

52,400.0000 MJD TDB 

Ai/ 

5.25 (12) x 10~ 7 Hz 

Ai// 

>8.7 x 10- 7 Hz 


14 1 (7) 

Ai/rf 

Av g + (~5 x 10- 9 ) Hz 


15.9 (6) days 

Ai> 

+2.18 (25) x 10-' 6 Hz s" 1 


52,443.13 (9) MJD TDB 

rms tuning residual 

44.9 ms 

Start observing epoch 

51,613 MJD 

End observing epoch 

52,900 MJD 


8 Numbers in parentheses represent 1 a uncertainties in the least significant 
digits quoted. 

b Lower limit given at 90% confidence. 
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(and correspondingly Ai/j) does not change the timescales 
significantly. 

The frequency evolution of IE 2259+586 just before and 
following file outburst as determined by our fit is shown in 
Figure 7. The preoutburst ephemeris is fully consistent with 
the ephemeris determined through earlier monitoring of this 
pulsar (Gavriil & Kaspi 2002). The glitch epoch (t g ) precedes 
the PCA observation containing the burst activity by 12.5 ± 
2.1 hr. Note that the glitch epoch precedes the observed burst 
activity whether the exponential growth tern is included in the 
fit or not (it is 4.7 a in the case with no growth term). The 
exponential growth term clearly improves the fit to the full 
data set; however, there is a systematic trend in the phase 
residual cluster just postglitch. This is discussed in detail 
below. Including the growth term significantly reduces the 
timescale of the exponential decay term to ~1 6 days. Finally, 
the long-term postglitch spin-down rate decreases signifi- 
cantly in magnitude (8 a) for both models. 

The parameters in Table 4 show that the glitch consisted of a 
total fractional frequency increase Aiw/V = (4.24 ± 0.1 1) x 
10 -6 , where Ai/^ is the maximum upward excursion in fre- 
quency relative to the preglitch ephemeris. Of this frequency 
jump, a fraction Q = (Au m - Av - Av g 4 - = 

0.185 ± 0.010 was recovered within 2 months following the 
glitch. 

We note that a decrease in the magnitude of the spin-down 
rate would be unique among all known pulsar and AXP glitches 
(Shemar & Lyne 1 996; Kaspi & Gavriil 2003). This is discussed 
more in § 4.2.5. Alternatively, some radio pulsars have shown a 
long-term exponential rise in frequency postglitch. However, 
our best fit to this model was excluded by the frequency data. 

The exponential rise term in this frequency model is clearly 
favored by the data; however, several important caveats must 
be stated. The exponential form was chosen because of its 
success in modeling radio pulsar glitches (Lyne et al. 1993). 
However, for the 19 days following the glitch, there is a sig- 
nificant deviation from this model that constitutes a large 
portion of the remaining x 2 i n the fit. During this 19 day 
interval, the best-fit model predicts rapid spin-up for the first 
~10 days, after which the frequency derivative again becomes 
negative. If we fit only the data from the 12 days following the 
glitch, we measure a frequency derivative of — (1.0 ± 0.3) x 


Fig. 7. — Top : Frequency evolution of IE 2259+586 around the time of die 
outburst for a model including an exponential rise and fall in frequency 
postglitch. See text for details. The line represents the best-fit model (see Table 4 
for model parameters). The circles denote frequency measurements for inde- 
pendent subsets of data. The effect of the glitch is obvious, as is the partial 
recovery. Middle : Frequency residuals of the independent frequency points 
minus die model. Bottom : Phase residuals with respect to the best-fit model. 
Closer inspection of the residual cluster just following the outburst epoch 
reveals a low-amplitude systematic trend. This is discussed further in the text 

10- 13 Hz s” 1 , whereas the model predicts an average fre- 
quency derivative of 5 x 10 -14 Hz s -1 . Summing the fre- 
quency derivative measurement error and the model error in 
quadrature, we find that the two values differ at the ~5 a 
level. Hence, there is no direct proof that there was significant 
spin-up during the first 12 days following the glitch. Note also 
that the frequency derivative measurement is at only the 3 a 
level, so some spin-up cannot be completely ruled out An- 
other caveat when we consider the exponential rise in fre- 
quency is that the pulse profile was undergoing large changes 
during this same time interval (see § 3.4), thus compromising 
our ability to phase align with our template pulse profile. In 
particular, if the pulse profile was shifting in phase in a smooth 
manner as it changed shape, then this shift in phase would 
manifest itself as an apparent change in frequency. It is not 
likely that the deviation from a pure exponential recovery can 
be attributed entirely to a systematic shift in the pulse profile, 
as that would require a large drift of ~0.35 cycles within the 
2 weeks following the glitch. Furthermore, even if we exclude 
from our fit the data in which the pulse profile changes were 
largest, the growth term still significantly improves the fit over 
the simple exponential recovery (Ax 2 = 136 for 2 fewer dof ). 
We conclude from our analysis that there is a significant de- 
viation from a simple exponential recovery; however, due to 
gaps in the data coverage (in particular from 19 to 60 days 
after the glitch) and the inherent pulse profile changes, we 
could not precisely identify the manner in which the frequency 
evolution deviated. 

3.2. Pulsed-Flux History 

Coincident with the burst activity on 2002 June 18 was a 
sudden increase in the pulsed flux from IE 2259+586 (Kaspi 
et al. 2003). The pulsed intensity of the AXP decreased through 
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Fig. 8. — Pulsed-flux history of IE 2259+586 (2-10 keV) from 2000 March 
through 2003 June as measured using die RXTE PCA. The sharp rise in die 
pulsed intensity coincides with the burst activity. 

the burst observation as did the burst rate, thereby confirming 
the AXP as the origin of the burst emission. Here we put the 
pulsed-flux enhancement in context with the long-term pulsed- 
flux history and track the recovery back toward the preoutburst 
level. 

Using the ephemeris determined in the last section, pulse 
profiles of the PCA data (2-10 keV) from 2000 March 
through 2003 June were constructed. The pulse profiles were 
generated for each pointing at times before and up to 1 year 
following the outburst. For the 2002 June 18 observation, the 
bursts were removed, and the data were split into ~200 s 
segments before folding. Between 1 and 19 days after the 
outburst, the data were grouped by spacecraft oibit (~3 ks). 
Once the pulse profiles were constructed, the pulsed intensity 
was measured by first decomposing each pulse profile in terms 
of its Fourier powers. The power in the first seven harmonics 
was summed to give the rms pulsed intensity (see eq. [1]). 
The pulsed-flux history of IE 2259+586 is shown in Figure 8. 
The large spike indicates the time of the outburst. Note that the 
pulsed flux has not yet returned to the preoutburst level. The 
timescale and functional form of the recovery are studied in 
more detail in § 3.6. 

3.3. Pulsed Fraction 

The PCA is a sensitive X-ray detector with a wide FOV and 
no imaging capabilities (for a fixed pointing). Because of the 
design of the instrument, it is not possible to reliably measure 
the pulsed fraction and/or spectrum of weak X-ray sources 
such as IE 2259+586 on account of uncertainties in the back- 
ground, particularly in the direction of the Galactic plane. 
However, during episodes when the AXP flux increases sub- 
stantially, one can measure accurately the pulsed fraction and 
spectrum (see § 3.5) with the PCA under certain assumptions 
(Lenters et al. 2003). 

The observed count rate in the PCA during pointings of IE 
2259+586 consists of several components, namely, the central 
source, instrumental background, near-Earth background due 
to precipitating particles, the “diffuse” cosmic background, 
the Galactic ridge, the bright SNR surrounding the AXP, and 


other dim sources in the FOV. To disentangle the spectrum of 
the central source from all the other contributions to the net 
count rate, the remaining components must somehow be 
modeled. The instrumental and near-Earth background (i.e., 
local background) can be subtracted using housekeeping data 
and models provided by the RXTE team. 11 Fortunately, the 
remaining components do not vary greatly on timescales of 
months to years. 

To estimate the count rate in the PCA of the remaining 
background components, we did the following. Using the last 
PCA observation (2002 May 4) taken before the outburst, we 
selected counts in the range 2-10 keV from all PCUs other 
than PCU0 and subtracted the background as estimated using 
the FTOOL PCABACKEST with the combined model version 
of 2003 March 30. Note that PCU0 was excluded because of 
the loss of the propane layer and the resulting large increase in 
the magnitude and variance of the background rate. Next, we 
measured the pulsed count rate (per PCU) during this obser- 
vation. Assuming that the 2-10 keV pulse fraction measured 
during the XMM-Newton observation from 2002 June 11 
(Table 2) was the same on 2002 May 4, we determined the 
background rate in the PCA. We use this “cosmic” back- 
ground rate for all observations during and following the 
outburst for which we measure the pulsed fraction. In using 
this count rate for our background, we assume that (1) the 
AXP pulse fraction was the same on 2002 May 4 as it was on 
2002 June 1 1 and (2) the cosmic background remained con- 
stant from 2002 May 4 until 2002 July 6 (the time of the last 
pulsed-fraction measurement reported here). 

Using this technique, we determined that the pulsed fraction 
of IE 2259+586 decreased relative to the preoutburst level, 
while the source was burst active (Kaspi et al. 2003). Here we 
have extended this analysis to later PCA observations through 
2002 July 6. As with the XMM-Newton data, we measured the 
pulsed fraction by subtracting the background and decom- 
posing the folded pulse profile in terms of its Fourier powers. 
The rms pulsed fraction was determined from the sum of the 
first seven harmonics by using the formalism described in 
§ 2.2. The RXTE PCA pulsed-fraction measurements are 
plotted, along with the XMM-Newton measurements (Table 2) 
in Figure 9. We found that the 2-10 keV pulsed fraction de- 
creased initially to a level of ~15% at which it remained for at 
least 1 day before rapidly returning to the preoutburst value 
within ~6 days of the outburst. 

3.4. Pulse Profile 

Coincident with the detected burst activity in the PCA from 
IE 2259+586 was a sudden change in the folded pulse profile 
(Kaspi et al. 2003). In the energy range 2-5 keV, the relative 
amplitudes of the two peaks were switched while the source 
was burst active. In this energy range, the pulse profile re- 
turned to near its preoutburst form within ~6 days. No 
changes like this have been seen in more than 6 years of 
monthly monitoring with the PCA (Gavriil & Kaspi 2002), 
although similar changes have been noted prior to the PCA 
monitoring (Iwasawa, Koyama, & Halpem 1992). Here we 
investigate further the pulse profile evolution over a longer 
time baseline and as a function of energy. 

The temporal evolution of the folded pulse profile of IE 
2259+586 from 2-10 keV is shown in Figure 10. From the 
pulsed-flux analysis, we know that the amplitudes of both 

11 See ftp://legacy.gsfc.nasa.gov/xte/calib_data/pca_bkgd. 
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Fig. 9. — Time evolution of the 2-10 keV pulsed fraction of IE 2259+586 
through the outburst of 2002 showing measurements made with the RXTE 
PCA 0 diamonds ) and XMM-Newton measurements ( squares ; Table 2). The 
horizontal dashed line marks the pulsed fraction from the XMM-Newton ob- 
servation 1 week prior to the outburst (Obs2). 

peaks increased during the 2002 June 18 observation relative 
to their preoutburst amplitudes. Furthermore, the pulsed flux 
decreased significantly during the burst observation ( RXTE 
oibits 1-3). This allows us to conclude that the amplitude of 
peak 1 is decaying more rapidly than peak 2 during the burst 
observation (as opposed to peak 2 growing relative to peak 1). 
The variability in the 18 days following the burst activity (see 
panels labeled “1.4 days” to “18 days”) shows erratic be- 
havior in the pulse morphology. It is not until several weeks 
after the outburst that the pulse profile closely resembles its 
preoutburst form within this energy range. Even at several 
months following the outburst the differences are significant 
Specifically, the bridge of emission connecting the two peaks 
is higher than it was preoutburst 

We quantified the change in pulse shape (2-10 keV) by 
decomposing the pulse profile in terms of its Fourier powers 
following equation (1). As is clear from Figure 11, the ratio of 
the power in the second harmonic to the first harmonic (i.e., 
fundamental frequency) is the dominant factor governing the 
pulse shape change. During quiescence, the second harmonic 
contains more than 80% of the total power. During the ob- 
servation containing the bursts, the power in the first harmonic 
increased such that there is actually as much or more power in 
this harmonic compared with the second harmonic. The ratio 
of the power in the first two harmonics recovers gradually 
over the next several months. 

Last, we investigated the dependence of the pulse shape on 
photon energy during the burst observation when the pulse 
shape was distorted the greatest. We find that there is a sig- 
nificant energy dependence of die pulse shape that becomes 
more prominent toward the end of the burst observation. 
Specifically, peak 2 becomes narrower at high energies, and 
the centroid of this peak lags in phase by ~0.02 cycles (0.15 s) 
from 2 to 15 keV. Both the centroid and phase of peak 1 are 
consistent with remaining unchanged versus energy. Similarly, 
the relative amplitudes of peak 1 and peak 2 are consistent 
with being constant in this energy range. 

3.5. Spectroscopy 

During the outburst detected on 2002 June 18, the X-ray 
spectrum of IE 2259+586 became much harder (Kaspi et al. 


2003). Initially, the blackbody temperature increased to 
~1.7 keV (from ~0.42 keV), and the photon index flattened 
to ~2.5 (from ~4.2). During the next 4 hr, the spectrum of 
IE 2259+586 evolved back toward its quiescent shape but did 
not recover fully. Here we analyze the RXTE PCA observations 
that took place over the next few weeks to track the recovery of 
the AXP spectrum. 

As with the pulsed fraction, measuring the X-ray spectrum 
of IE 2259+586 with the RXTE PCA is not straightforward, so 
we used a similar technique in subtracting the background as 
that employed to measure the pulse fraction (see § 3.3). As 
before, the last PCA observation carried out on 2002 May 4, 
46 days before the outburst, was used to estimate the cosmic 
background in the PCA FOV. First, the local background 
during this observation was estimated using the FTOOL 
PCABACKEST and the combined model of 2003 March 30. 
The observed Standard 2 spectrum 12 was grouped such that 
there were at least 25 total counts per bin (source plus back- 
ground) and then fitted (less the local background) using 
XSPEC version 1 1 .2.0 to a multicomponent model. The model 
consisted of the standard AXP spectrum plus another black- 
body and a Gaussian line. The second blackbody component 
and Gaussian line were each modified by interstellar absorp- 
tion with a fixed column density (Nn = 2 x 10 22 cm -2 ). The 
AXP spectral parameters were frozen to the values obtained 
from the XMM-Newton observation of 2002 June 11 (see 
Table 3), 37 days following this observation and 7 days prior 
to the outburst. Using this model, we obtained a good fit to the 
PCA data between 2.5 and 20.0 keV (x£ = 0.71 for 35 dof). 
The free components in this fit (i.e., blackbody and Gaussian 
line), therefore, define the spectrum of the remaining back- 
ground in tiie PCA for this particular pointing at this particular 
epoch. This model was used to define the cosmic background 
in the PCA detector during and directly following the burst 
activity, assuming that these components (e.g., Galactic ridge, 
SNR, etc.) remain constant between 2002 May 4 and July 6. 

A Standard 2 spectrum was accumulated for the 2002 June 
18 observation when the bursts were detected. Note that the 
bursts themselves were eliminated from the data used to create 
the spectrum. No single-component model provided an ade- 
quate fit to the data. For example, a fit to a power-law model 
yielded a statistically unacceptable x 2 of 384.7 for 37 dof. 
Conversely, the standard blackbody plus power law model 
resulted in a very good fit (x 2 = 35.8 for 35 dof). We measure 
a significantly harder spectrum than ever before seen in this 
AXP with kT = 1 .23 ± 0.4 and T = 2.83 ± 0.07. 

Next, we split the 2002 June 18 into seven separate seg- 
ments to search for time evolution of the spectrum. Spectra 
were also accumulated for each PCA observation from 2002 
June 20 through July 6. We fitted the blackbody plus power 
law model to each of these PCA observations (2.5-20.0 keV). 
The RXTE fit results are shown in Figure 12 along with the 
XMM-Newton measurements given in Table 3. 

We tested the robustness of our technique by varying the 
assumptions that the spectrum of IE 2259+586 remains con- 
stant during the 44 days leading up to the burst activity. Speci- 
fically, when fitting for the PCA cosmic background in the 
2002 May 4 observation, we tested a range of values for the 
spectral parameters (kT and T) of IE 2259+586 within their 
respective historical ranges (see § 2.3). Obviously, this did 
alter the cosmic background fit parameters but resulted in 

12 As with the pulsed-fraction analysis, all data from PCU0 were excluded 
from the spectral analysis. 
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Fig. 10. — Time evolution of the 2-10 keV pulse profile of IE 2259+586 through the outburst of 2002. Profiles have arbitrary flux units, and two pulse cycles are 
shown for clarity. Time increases from left to right and top to bottom. 


only insignificant changes (<1%) to the measured spectral 
parameters of the AXP for the 2002 June 18 observation. For 
the final PCA observation fit in this analysis for which the 
AXP was dimmest (2002 July 6), the relative change in these 
parameters is larger (~5%) but is still less than 1 a. 

The unabsorbed flux decayed rapidly within the first day of 
the outburst. The functional form of the decay depends criti- 
cally on the reference epoch chosen. Both power-law and ex- 
ponential decay models yield acceptable fits to the data. At 20 
days following the outburst, the AXP is still a factor of ~2 
brighter than its nominal (quiescent) flux level. The flux decay 
is covered in greater detail in the next section. The remaining 
spectral parameters shown here recover to within 25% of their 
preburst levels within the first ~l-3 days. By using the 
pulsed-flux history as an indicator for the recovery timescale, a 
full recovery of the spectral parameters would not be expected 
until ~1 yr later. However, variation in these spectral parameters 
at the ~25% level has been seen outside burst activity (see 
§ 2.3); thus, these parameters may not decay further. 

It is interesting to note that the X-ray spectrum of IE 
2259+586 was harder when the spin-down rate was higher 


during the first month following the glitch. This behavior in 
IE 2259+586 is qualitatively consistent with the spectral 
hardness versus spin-down rate correlation found by Marsden 
& White (2001), who considered the AXPs and SGRs as a 
whole. This suggests that individual members of this class that 
have shown significant variability in spin-down rate may also 
show correlated spectral variations. 

We note that the later RXTE PCA spectral measurements 
(>1 day) agree reasonably well with the XMM-Newton mea- 
surements contained within this interval; however, they are 
systematically offset from one another. This could be due to 
systematic effects in our background subtraction method, a 
deviation in the power-law spectral model at high energies 
at which the two instrumental responses do not overlap 
(12-20 keV), and/or a systematic offset in the cross- 
calibration between the XMM-Newton pn and the RXTE PCA. 

3.6. Flux Decay and Energetics 

To better quantify the flux decay of IE 2259+586 fol- 
lowing the outburst of 2002 June, we combined the XMM- 
Newton flux measurements with our RXTE PCA pulsed-flux 
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measurements, which span a much broader temporal baseline. 
The advantage of using the PCA pulsed-flux measurements as 
opposed to the phase-averaged flux measurements is that the 
systematic errors present in the background subtraction are not 
a concern. 

The pulsed-flux measurements were converted to unab- 
sorbed phase-averaged fluxes in the following way. First, we 
determined a conversion factor between the two by calculating 
the ratio of the preoutburst unabsorbed flux measured with 
XMM-Newton to the average PCA pulsed flux for the year 
leading up to the outburst. Assuming that the pulsed fraction 
and spectrum do not change, this factor can be multiplied by 
subsequent PCA pulsed-flux values to estimate the unab- 
sorbed phase-averaged flux at those times. However, we know 
that both the pulsed fraction and the energy spectrum changed 
during this outburst. We corrected for the brief period during 
which the pulsed fraction decreased by multiplying by an 
additional factor of the ratio of the nominal pulsed fraction to 
the observed pulsed fraction at those times. This factor was 
applied only to the PCA pulsed-flux measurements within the 
first 2 days following the burst activity during which the 
fraction dropped from 23% to ~1 5%. Computing this con- 
version factor for a broad energy range (2-10 keV) reduces 


the effects of spectral changes. In fact, when we applied our 
pulsed flux-to-unabsorbed phase-averaged flux conversion 
factor to the data for which we have independent unabsorbed 
flux measurements at early times in the outburst (<20 days 
after burst activity), the agreement between the two mea- 
surements was found to be quite good (Fig. 13). Since the 
spectral differences are greatest at these times in the outburst, 
we feel that this technique is a robust one for constructing a 
long-term light curve for the source. 

In this way, we have taken the pulsed-flux measurements 
presented in Figure 8 and converted them to phase-averaged 
unabsorbed flux values. These data are plotted relative to the 
glitch epoch (Table 4) in Figure 13 in addition to the three 
XMM-Newton flux measurements. Clearly, the flux decay is 
not well described by a single-component model (e.g., expo- 
nential or power law). The temporal decay of the flux during 
the PCA observations containing the burst activity (<1 day) is 
much more rapid than the decay during the year following the 
burst activity. We split the data into two segments (<1 day 
and >1 day after the glitch) and fitted each independently to a 
power-law model (F ex f*). The measured temporal decay 
indices forthe two segments are —4.8 ± 0.5 and —0.22 ± 0.01, 
respectively. 
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Fig. 12. — Spectral evolution of IE 2259+586 through and following the 
outburst of 2002. Top to bottom: Unabsorbed flux (2-10 keV), the blackbody 
temperature (kT), the photon index, the blackbody radius, and the ratio of 
power-law (2-10 keV) to bolometric blackbody flux. A distance of 3 kpc was 
assumed (Kothes, Uyaniker, & Aylin 2002) to compute the blackbody radius. 
Horizontal dashed lines denote the values of each parameter during theXMM- 
Newton observation 1 week prior to the outburst. 


The RXTE All-Sky Monitor (ASM) observed IE 2259+586 
at 0350 and 1443 UT on 2002 June 18, before the observed 
burst activity and following the measured glitch epoch, al- 
though the first observation is within the errors of the measured 
glitch epoch. If we extrapolate the slower flux decay model 
(a = —0.22) back to the time of the ASM observations, we 
find that the expected flux values (14 and 8 x 10 -11 ergs cm -2 
s _1 , respectively) are comfortably below the 99% confidence 
upper limit of 1 x 10~ 9 ergs cm -2 s' 1 (Kaspi et al. 2003); thus, 
the ASM limits are not constraining for this component. For 
the steeper component (a = — 4.8) containing the burst ac- 
tivity, we find that the expected flux at the time of the first ASM 
observation is several orders of magnitude above the upper 
limit. The second ASM observation (52,443.61 1 MJD TDB) is 
unconstraining. This suggests that the onset of burst activity 
(associated with the first decay component) must have come 
after the earlier ASM observation (52,443.158 MJD TDB) or 
the flux enhancement deviated from this steep decay between 
the glitch epoch and the RXTE PCA observations. One possi- 
bility is that the onset of this burst activity was delayed with 
respect to the glitch epoch. In fact, a much later onset time for 
the burst activity is inferred from the time evolution of the burst 
recurrence frequency (Gavriil et al. 2003). 

The absence of a fiducial point for the initial, rapid flux 
decay associated with the burst activity does not allow us to 


accurately measure the temporal decay index of this compo- 
nent. Choosing a reference epoch just before the time of the 
first bursts detected with the PCA yields a decay index less 
than unity (in magnitude). Thus, we can only constrain the 
index of the flux decay within the time span containing the 
burst activity to be less than 6.6 (in magnitude). Note that as 
the reference epoch approaches the start of the burst obser- 
vation, the flux decay becomes steeper than a power law in 
form (e.g., exponential). This is discussed further in § 4.1.1. 
Unlike the early component, the temporal index of the more 
gradual flux decay in the months following the glitch 
(a = —0.22) is insensitive to varying the reference epoch 
between the time of the glitch and the beginning of the ob- 
served burst activity. 

We next used these power-law fits for estimating the total 
energy released during this outburst. We integrated each 
model, less the quiescent flux level (Fig. 13), only over the 
time ranges during which we have flux measurements (i.e., 
from the start time of the burst observation onward). In spite 
of the large range of allowed temporal decay indices, the 
energy we measure in the early flux decay component is well 
determined since we integrate the model only over the interval 
for which we have observations (i.e., we do not extrapolate the 
model back toward the glitch epoch). Assuming a distance of 
3 kpc to IE 2259+586 (Kothes, Uyaniker, & Aylin 2002), we 
measure an energy release (2-10 keV) of 2.7 x 10 39 and 
2.1 x 10 41 ergs for the fast and slow-decay intervals, respec- 
tively. In terms of the overall energy budget, the energy re- 
leased in the bursts themselves (6 x 10 37 ergs over 2-60 keV; 
Gavriil et al. 2003) is insignificant in comparison with the 
excess persistent emission in X-rays released during the year 
following the outburst. Moreover, the excess persistent flux 
emitted during the interval containing the burst activity is 
insignificant in comparison with the total energy released 
during the year-long flux decay. 



Fig. 13. — Time evolution of the unabsorbed flux from IE 2259+586 fol- 
lowing the 2002 June outburst. The glitch epoch (Table 4) is used as the 
reference time for this plot. Diamonds denote inferred unabsorbed flux values 
calculated from RXTE PCA pulsed-flux measurements. Asterisks and squares 
mark independent phase-averaged unabsorbed flux values from RXTE and 
XMM -Newton, respectively. The dotted line denotes the flux level measured 
using XMM-Newton 1 week prior to die glitch. The dashed line is a power-law 
fit to the PCA flux measurements during the observations containing the burst 
activity (<1 day). The dot-dashed line marks the power-law fit to all data more 
than 1 day following the glitch. See text for further details. 
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4. DISCUSSION 

Virtually all measurable X-ray properties of IE 2259+586 
changed suddenly and dramatically during the 2002 June 
outburst Continued observations with RXTE and XMM- 
Newton have allowed us to track the recovery of several source 
parameters shown to change during this outburst (Kaspi et al. 
2003) and identify additional parameters that were similarly 
affected. Many of the observed variations resemble phenom- 
ena seen in other classes of neutron stars, namely, SGRs and 
radio pulsars. Here we compare the AXP outburst properties 
with similar phenomena seen in those source classes in the 
hope of identifying similarities and differences that can help 
elucidate the physical properties of these different manifes- 
tations of young neutron stars. Specifically, we consider (1) 
the radiative properties of the persistent and pulsed emission 
during and following die outburst and compare these with 
those seen in SGRs, since that is the only other source class 
for which such outbursts have been seen, and (2) the rotational 
behavior of the pulsar and compare it with behavior seen in 
radio pulsars, as well as in SGRs, and in another AXP, 1RXS 
J 1708-4009. Finally, we point out that the detection of low- 
intensity outbursts in AXPs has important implications for our 
estimates of the number of active magnetar candidates in our 
Galaxy. The IE 2259+586 burst properties and their relation 
to SGR burst properties are considered separately (Gavriil 
et al. 2003) 

4.1. Transient X-Ray Emission and Pulse Properties 

We have shown that there are two components to the flux 
decay in IE 2259+586 during die 2002 June outburst There is 
a rapid decay of the flux during die observations containing the 
burst activity, and then there is a more gradual flux decay seen 
in the year following the glitch. During die initial flux decay, 
the spectrum was considerably harder than at all other times, 
while the spectral hardness at times greater than 1 day after the 
glitch is consistent with preglitch spectral measurements. The 
spectral differences in the two decay components (total energy, 
blackbody radius, etc.) point toward separate physical mech- 
anisms for the two flux enhancements (see Figs. 9 and 12). 

Coincident with the glitch and burst activity was a dramatic 
change in the pulse profile of IE 2259+586. The majority of 
the observed shape change recovered within ~6 days, but there 
was still some residual change that slowly decayed over the 
months following the outburst (§ 3.4). The 2-10 keV pulsed 
fraction dropped to ~15% during the observed burst activity 
and quickly recovered to the preoutburst level of 23.4% within 
~6 days (§ 3.3). Here we compare the observed properties of 
the flux enhancement and pulse properties in IE 2259+586 
with qualitatively similar behavior detected in SGRs and 
briefly discuss how the changes seen in IE 2259+586 can be 
accommodated within the magnetar model. 

4 . 1 . 1 . Comparison with SGR Outbursts: X-Ray Flux and Spectrum 

The richest SGR database with which to make an empirical 
comparison with the IE 2259+586 outburst comes from the 
most active SGR during the last 6 years, SGR 1900+14. This 
SGR has been observed on 14 separate occasions since 1997 
by imaging X-ray telescopes and more than 100 times with 
RXTE . Within this time span, SGR 1900+14 entered several 
burst-active episodes; the most notable of which was the 
outburst that began on 1998 August 27 with a giant flare 
having a total energy MO 44 ergs (e.g., Feroci et al. 2001). 
Coincident with this giant flare was a large increase in the 


persistent and pulsed flux from the source (e.g.. Woods et al. 
2001), in addition to a dramatic change in the pulse profile and 
a timing anomaly (see § 4.2.3). 

A spectral analysis of the RXTE PCA observation of SGR 
1900+14 1 day following the giant flare shows that the black- 
body temperature was higher than the nominal temperature 
(~0.5 keV) at 0.94 keV (Woods 2003). Two other high- 
fluence bursts from SGR 1900+14 have extended X-ray tails 
or afterglows that show enhanced thermal emission at times 
reaching up to 4 keV (Ibrahim et al. 2001; Lenters et al. 2003). 
The thermal component of IE 2259+586 shows a similar 
brightening and temporal decay. Here the temperature rose 
from 0.4 up to 1.7 keV at the onset of the outburst before 
quickly decaying to 0.5 keV within the first few days. Ac- 
companying the temperature increase in IE 2259+586 was a 
significant hardening of the photon index. This is different 
from what has been seen in SGR 1900+14, where there has 
been either no change in photon index after some bursts 
(Lenters et al. 2003) or even a slight softening of the non- 
thermal component of the spectrum (Woods 2003). 

An analysis of four separate bursts or flares from SGR 
1900+14 and their associated afterglows (Lenters et al. 2003) 
shows that the emitted energy within the afterglow (2-10 keV) 
corresponds to roughly 2% of flue burst energy (25-100 keV). 
Because of the bandpass of the PCA, we have been able to 
measure the burst energy only within the 2-60 keV range 
(Gavriil et al. 2003). However, we can estimate the burst en- 
ergy (25-100 keV) during the 2002 June outburst by taking 
the measured count fluence with the PCA (Gavriil et al. 2003) 
and multiplying by the counts-to-energy conversion factor 
determined for SGR 1900+14 (Gogu$ et al. 1999). This con- 
version factor was used because the mean burst spectral 
hardness in the PCA was similar for IE 2259+586 (Gavriil 
et al. 2003) and SGR 1900+14 (Gogu$ et al. 2001). The total 
energy was further modified by a factor-of-2 increase to 
roughly compensate for gaps in the data due to Earth occul- 
tation during the PCA pointing. Assuming that the energy 
spectra of the IE 2259+586 bursts above 60 keV is similar to 
SGR 1900+14, we estimated the burst energy released by IE 
2259+586 to be ~1 x 10 38 ergs (25-100 keV). The brighten- 
ing of IE 2259+586 during the 2002 June outburst contained 
two components: a rapid decay during the initial burst ob- 
servation followed by a more gradual flux decay in the year 
following the glitch. The energy released from IE 2259+586 
within these two components is ~30 and ~2000 times greater 
than the burst energy (Table 5), in stark contrast with the value 
of 0.02 found for SGR 1900+14. In feet, for a 2% after- 
glow-to-burst energy ratio, we can rule out the presence of an 
intermediate burst from IE 2259+586 preceding the 2002 June 
18 PCA observation. Assuming that only the hard component 
of the decay constituted the burst afterglow, the energy of the 
hypothetical burst powering this afterglow must be 1.5 x 
10 41 ergs (equivalent fluence of 1.4 x 10 -4 ergs cm -2 ). The 
Konus detector aboard the Wind spacecraft located at the LI 
point between Earth and the Sun maintained continuous 
coverage of IE 2259+586 from the time of the glitch through 
the first PCA observation and did not detect any emission 
from the AXP (S. Golenetskii 2003, private communication). 
The upper limits on the burst fluence (17-70 keV) derived 
from the Konus data are 1.7 x 10' 7 ergs cm" 2 for bursts of 
duration less than 3 s and 1.7 x 10 _7 (At/3 s) 1 / 2 ergs cm -2 , 
where At is the burst duration, for bursts lasting between 3 
and ~1000 s. Clearly, the putative burst with fluence 1.4 x 
1 0 -4 ergs cm -2 can be excluded for durations less than ~1 000 s. 
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TABLE 5 

SGR and AXP Burst/Afterglow Energetics and Temporal Decay Indices 


Parameter 

SGR 1900+14* 

SGR 1627-41 

IE 2259+586 b 

IE 2259+586 c 

Burst energy (ergs; 25-100 keV) 

1 xlO 44 

4 x 10 42 

1 x 10 38 

1 x 10 38 

Tail energy (ergs; 2-10 keV) 

2x 10 42 

3 x 10 42 

3x 10 39 

2x 10 41 

Decay index 

-0.71 

-0.47 

>-6.6 

-0.22 


Note. — The following distances are assumed for conversions to energy: 15 kpc for SGR 1900+14 (Viba et al. 2000), 11 kpc for 
SGR 1627-41 (Corbel et al. 1999), and 3.0 kpc for IE 2259+586 (Kothes et al. 2002). 

a The values given here correspond only to those measured for die 1998 August 27 flare and its associated afterglow. 
b The tail energy given here is derived from the excess persistent emission observed during the burst observations only. 
Because of the ambiguities in the epoch determination, an accurate decay index cannot be measured. See § 3.6 for details. 

c The tail energy given here is derived from the excess persistent emission observed during the long timescale decay. See § 3.6 
for details. 


The temporal decay of the flux from SGR 1900+14 fol- 
lowing the August 27 flare follows a power law in time with 
an exponent of -0.71 during the 40 days following the flare 
(Woods et al. 2001). For a power-law fit to the flux evolution 
during the burst activity, we can constrain the decay only to 
have an exponent less than 6.6 (in magnitude) because of the 
uncertainty in the reference epoch (see § 3,6). Thus, this decay 
component is not inconsistent with a —0.7 power law. The 
more gradual flux decay in the months following the glitch 
obeyed a power law in time that scaled as f 0 - 22 * 00 ^ signif- 
icantly flatter than the decay following the August 27 flare 
or any other flux decay following bright bursts from SGR 
1900+14 (Lenters et al. 2003; Feroci et al. 2003). 

We can also compare the energetics and flux decay of IE 
2259+586 with those seen in SGR 1627-41. This SGR has 
shown one outburst in 1998 during which 98% of the emitted 
burst energy was concentrated into a narrow 3 week window 
(Woods et al. 1999a). The flux from the source decayed over 
the next 3 years approximately as a power law in time with an 
exponent of -0.47 (Kouveliotou et al. 2003). The last two 
Chandra observations of this SGR show that the flux has 
leveled out (at least temporarily). Since there are no pre- 
outburst flux measurements to establish the “quiescent” flux 
level of this SGR, we take the two latest flux measurements as 
the quiescent flux level. Under these assumptions, the burst 
tail or afterglow energy during the 4 years following the 1998 
outburst is comparable to the energy released in the bursts 
themselves (Kouveliotou et al. 2003). 

The burst and persistent emission energetics and temporal 
decay indices for all three sources, IE 2259+586, SGR 1900+ 
14, and SGR 1627-41, are listed in Table 5. Clearly, there are 
large differences in both the decay index and the ratio of burst 
to persistent emission energetics among the three sources. One 
possibility is that the response of the persistent X-ray flux to 
(or when accompanied by) burst activity in SGRs and AXPs 
varies within the group. Alternatively (or in addition), there 
may be two components to the SGR flux decay similar to the 
AXP flux evolution, and the values listed in Table 5 for the 
SGRs reflect a mixture of the two components that contributes 
to the quantitative differences with the AXP values. In fact, 
the flux from SGR 1900+14 in the months following the end 
of the 40 day afterglow was enhanced relative to the pre- 
outburst level (Woods et al. 2001). This enhancement could be 
due to the persistent low-level burst activity observed during 
this time interval or perhaps due to a slower flux decay 
component, analogous to the slow flux decay seen in IE 
2259+586. Unfortunately, there are no reported spectral 
measurements of SGR 1900+14 during this epoch to help 
distinguish between the two possibilities. For SGR 1627-41, 


the initial follow-up observation was not performed until ~50 
days after the primary outburst (Woods et al. 1999a); there- 
fore, any short-lived transient flux decay would have been 
missed. Thus, it appears that two flux components may pos- 
sibly be present in all SGR and AXP outbursts. More rigorous 
follow-up spectral measurements after future outbursts are 
needed to show whether or not this behavior is ubiquitous. 

4.1.2. Comparison with SGR Outbursts: Pulse Properties 

SGR 1900+14 has shown one clear instance of a correlated 
change in the pulse profile and pulsed fraction following the 
burst of 1998 August 29 (Lenters et al. 2002). Here the pulsed 
fraction increased to a maximum value of ~20% at ~200 s 
following the burst then rapidly decayed back to the preburst 
value of ~12% within 10 4 s. During the first ~100 s, the pulse 
shape showed large changes relative to the preburst profile, 
similar to the pulsed-fraction recovery, then recovered fully 
(within the errors) during the next 10 4 s. In IE 2259+586, we 
observe comparable changes in both the pulse shape and pulsed 
fraction at early times (£6 days) during the outburst. However, 
the timescale is significantly longer for IE 2259+586, and the 
pulsed fraction decreased rather than increased. 

The X-ray pulse profile of SGR 1900+14 changed from a 
complex multipeaked shape before the giant flare of 1998 
August 27 to a simple nearly sinusoidal profile following the 
flare (Woods et al. 2001). This change has persisted for years 
following the outburst and has yet to recover (Gogii§ et al. 
2002). The constraints on any change in pulsed fraction in SGR 
1900+14 following this flare are weak, as the first reported 
pulsed-fraction measurement was made at 19 days following 
the flare, consistent with the preflare value (Woods et al. 2001). 
In IE 2259+586, there was a significant change in the pulse 
profile in which the power in the fundamental frequency in- 
creased relative to the higher harmonics, analogous to what 
was observed in SGR 1900+14. The difference in pulse shape 
from before to during or following the burst activity was not 
as profound as the change seen in SGR 1900+14; however, 
the burst activity in IE 2259+586 was not nearly as ener- 
getic, either. Note also that the change in pulse profile of IE 
2259+586 was transient, as the pulse shape at 1 year following 
the outburst is very similar to the preoutburst pulse profile. 

4.1.3. Physical Interpretation 

In spite of the quantitative differences mentioned above, 
many of the features of the outburst in IE 2259+586 are 
qualitatively very similar to those seen in SGR outbursts, 
specifically, the flux enhancement, spectral change, and cor- 
related change in pulse properties. The similarities outlined 
above further solidify the connection between SGRs and 
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AXPs. Combined with the similar burst characteristics (Gavriil 
et al. 2003), this outburst in IE 2259+586 shows beyond any 
reasonable doubt that SGRs and AXPs are of the same nature, 
as predicted uniquely by the magnetar model. 

During the rapid initial flux decay, the blackbody radius 
was smaller (~1 km) than at all other times, and die tem- 
perature was higher (0.8- 1.7 keV). The thermal component of 
the spectrum suggests the existence of a hot spot that either is 
cooling through its surface or is being heated by energetic 
particles accelerated in the magnetosphere. A localized hot 
spot will clearly result in a change in the emitted radiation 
pattern (i.e., pulse profile) even if the pulse shape is strongly 
modified by scattering in the magnetosphere (Thompson, 
Lyutikov, & Kulkami 2002). The reduction in pulsed fraction 
suggests that the heated region is offset in angle from the 
locations on (or above) the stellar surface that give rise to the 
two pulse maxima. 

These observations do not allow us to distinguish con- 
vincingly between magnetospheric emission and passive cool- 
ing of an impulsively heated crust as the underlying source of 
the transient X-ray emission. The rapid flux decay could have 
a very different time scaling from what is observed following 
intermediate-energy bursts in SGR 1900+14, and any change 
in the magnetic field of IE 2259+586 may have occurred too 
gradually to generate a bright X-ray outburst A large current 
density will be excited in the magnetosphere above regions of 
strong crustal shear, but the mechanism by which this current is 
damped depends on how rapidly it is excited (Thompson & 
Duncan 2001). The rapid flux decay as measured is consistent 
with a gradually diminishing creep within a small region of 
the crust over a period of M0 5 s. 

It is also worth considering whether this part of the X-ray 
transient is the direct aftermath of a more energetic burst The 
short-lived afterglows detected after SGR bursts of interme- 
diate energy have a simple explanation as the cooling of a 
pair-rich surface layer heated by a high-energy burst (Ibrahim 
et al. 2001) and as such are valuable probes of the burst 
mechanism. As estimated above, die minimum energy of such 
aburstis ~3.0 x 10 41 (cag/10 -2 )' 1 ergs, where is afterglow 
energy divided by the burst energy. Konus data (S. Golenetskii 
2003, private communication) sets an upper bound of ~2 x 
10 38 ergs for bursts with duration less than 3 s and ~2 x 
10 38 (At/3 s) 1/2 ergs for bursts with duration 3 < At < 1000 s. 
Therefore, a much higher afterglow efficiency is required for 
IE 2259+586 than is inferred for the intermediate and giant 
SGR 1900+14 flares. 

Is there any evidence for such extended energy release 
during SGR outbursts? SGR 1900+14 has, in fact, been ob- 
served to radiate the same energy over timescales ranging 
from ~1 00 s in the declining tail of the August 27 flare (Feroci 
et al. 2001) down to ~2 s in the August 29 burst that followed 
it. Thus, the putative burst from IE 2259+586 would neces- 
sarily have a much longer duration than any (of an equivalent 
energy) yet observed from an SGR. Although the magneto- 
spheric opacity scales as B~ 2 and is expected to be lower in IE 
2259+586 than in the actively bursting SGRs, the luminosity 
of a passively cooling trapped fireball would almost certainly 
exceed the bound derived from Konus (e.g., for B ~ 2?qed)- 

We now comment on the relative importance of crustal 
heating and enhanced magnetospheric emission for the slow- 
decay component of the outburst of IE 2259+586. Bulk 
heating of the crust of a magnetar can power an excess 
heat flux from its surface for a year or longer and has been 


proposed as the explanation for the quasi-power-law flux de- 
cay seen in SGR 1900+14 (Lyubarsky, Eichler, & Thompson 
2002) and SGR 1627—41 (Kouveliotou et al. 2003). In each 
case, an initial energy deposition of 10 44 ergs was assumed. 
For SGR 1 900+14, there was in fact a giant flare preceding the 
flux enhancement that emitted MO 44 ergs in soft gamma-rays. 
For SGR 1627—41, the burst energy output was more than 1 
order of magnitude less, so it was suggested that the process 
that generates the burst energy was much less efficient in this 
source (Kouveliotou et al. 2003). As discussed previously, the 
relative energy released in burst and excess persistent emis- 
sion remains very uncertain in the case of IE 2259+586. Still, 
the excess persistent energy output in the slow-decay com- 
ponent (2 x 10 41 ergs) was only a factor of 10 lower than the 
output from the two SGRs, yet there was not a single burst 
detected with an isotropic luminosity above ~1 x 10 39 ergs s -1 
during this outburst. Thus, if we ascribe file slow flux decay 
seen in these three sources to deep heating of the crust of a 
magnetar, then it follows that IE 2259+586 is almost certainly 
less efficient in producing burst emission than either of the 
other two SGRs. Note also that the similarity in the timescales 
for the relaxation of the torque and the X-ray pulse profile 
must be coincidental if the slow-decay component of the 
X-ray emission is powered by crustal cooling. 

A small twist of the stellar crust, as is hypothesized to 
explain the glitch characteristics (§ 4.2.4), will impart a twist 
in the field lines that are anchored to this patch of crust 
This drives a current along these twisted field lines, which 
ultimately produces X-ray emission whose luminosity is pro- 
portional to the twist angle (Thompson et al. 2000, 2002). 
However, following a twisting motion of the crust and ex- 
ternal magnetic field through an angle 9 , at most a fraction 
~0.03(«/0.01Xflmax/10- 3 P(fi dip /10 14 G) 2 of the energy re- 
leased is stored in the external nonpotential magnetic field. 
(Here #max is the maximum strain that the crust can sustain 
before yielding.) The majority of the energy is dissipated 
through the deformation of the crust itself, because the 
shear modulus at the base of the crust corresponds to a much 
stronger magnetic field, B = (Any) 1 ! 2 ~ 6 x 10 15 G. The 
amplitude of the required global twist (about ~10 -2 rad) 
will only slightly change the optical depth to resonant cyclo- 
tron scattering or change the external torque through a flaring 
out of the external dipole field (Thompson et al. 2002). Hence, 
this cannot explain the large shape changes seen in the 
pulse profile within the first week following the glitch. It 
is possible, however, that the more subtle pulse profile 
changes seen more than 1 week postglitch can be explained as 
a slight change in optical depth due to a twisting of (likely 
extended) field lines. The torque is more sensitive to a current 
localized on the most extended field lines, as is discussed in 
§ 4.2.5. 

4.2. Rotational Evolution 

The timing data clearly show that a large spin-up glitch 
occurred in IE 2259+586 coincident with or perhaps shortly 
before the onset of burst activity. A portion of the frequency 
jump (~20%) recovered in a quasi-exponential manner within 
the next ~60 days. The long-term timing in the year following 
the glitch shows a significant reduction in the spin-down rate 
of this pulsar. Below, we draw comparisons between the 
characteristics of the IE 2259+586 glitch and those of other 
neutron stars and offer possible explanations for the observed 
behavior in IE 2259+586. 
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4.2.1. Comparison with Radio Pulsars: Overall Energetics 

A conservative bound on the change in rotational energy is 
obtained by treating the star as a rotating solid body. For a star 
having moment of inertia 7 = 10 45 g cm 2 , one infers a change 
in the rotational energy of 3 x 1 0 39 ergs, some 2 orders of 
magnitude smaller than the integrated excess X-ray emission. 
If the glitch involves an exchange of angular momentum be- 
tween a superfluid with moment of inertia 7 S f , and the crust of 
the star with moment of inertia 7 C , then the release of rotational 
energy is A E g ~ (27r) 2 7 c Ai/(i/ s f — v c \ where u s f — v c is the 
equilibrium lag between the rotation frequency of the super- 
fluid and crust. Evaluating this lag by using the observed 
~15 — d postglitch relaxation time, one finds a value for A E g 
that is smaller by a factor of 10 -3 (7 S f/0.027) _1 . Thus, the 
rotation of the star is manifestly not the source of energy for 
the transient X-ray emission. 

The relaxation of elastic strains in the stellar crust is another 
possible source of energy. The magnitude of these strains 
could easily exceed the equilibrium rotational bulge of 
the star, (/feqtr - tfpo. e )/*Ns = ftXs/GM NS = 10- 8 (/>/7 s)~ 2 . 
However, the energy stored in the strained crust is ~3 x 
10 39 (i/;/10~ 4 ) 2 ergs. The magnitude of the strain must, by way 
of comparison, be smaller than 5 x 10 -5 in the precessing 
pulsar PSR B 1828— 11 (Cutler, Ushomirsky, & Link 2003), 
whose spin is an order of magnitude faster than that of IE 
2259+586. We conclude that the observed X-ray transient 
could be powered by the release of elastic energy only if 
departure from spherical shape in the crust were comparable 
to that of a neutron star with a ~15 ms spin period. 

The remaining possibility is that the transient X-ray emis- 
sion is powered by magnetic field decay. The energy carried 
by the (exterior) magnetic field of IE 2259+586 is comfort- 
ably larger: 8 x 1 0 44 (Spoi e / 1 0 1 4 G) 2 (/? N s / 1 0 km) 3 ergs in the 
simplest case of a centered dipole. This is clearly a lower 
bound to the total magnetic energy of the star, as it does not 
account for the internal (e.g., toroidal) field. The surface field 
(and magnetic energy) would also be significantly stronger, by 
a factor of ~(A7?/2/? N s)“ 3 , if the dipole moment were offset 
from the center of the star (with a separation A R between the 
magnetic poles). 

4.2.2. Comparison with the AXP 1RXS Jl 708—4009 

It is interesting to compare the IE 2259+586 glitch with 
those reported for a different AXP, 1RXS 1708-4009. This 
AXP has shown two glitches, one in 1999 (Kaspi, Lackey, & 
Chakrabarty 2000) and one in 2001 (Kaspi & Gavriil 2003; 
Dall’Osso et al. 2003). The first glitch was similar to those 
seen in the Vela and other radio pulsars; it showed a step 
Au/u = 6 x 10“ 7 with Q = 0 and Ai>/z> ~ 1%, such that the 
magnitude of the spin-down increased. The 2001 glitch, 
however, was markedly different and in some ways resembled 
the glitch seen in IE 2259+586. Specifically, it showed 
Av/v = 1.4 x 10~ 7 , which completely recovered exponentially 
on a timescale of 50 days (i.e., Q « 1). Whether it suffered a 
step in v depends on the interpretation of an apparent v 
detected between glitches; this is discussed in more detail by 
Kaspi & Gavriil (2003) and Dall’Osso et al. (2003). As argued 
by Kaspi & Gavriil (2003), this glitch may have indicated that 
1RXS 1708—4009 underwent bursting activity some time in 
the few week long interval between monitoring observations 
that straddled the glitch epoch but went unobserved. However, 
no flux or pulse profile changes comparable to those seen in 
IE 2259+586 still weeks after its outburst were seen for 1RXS 


1708—4009. Thus, any radiative changes would have had to 
decay on shorter timescales than seen in IE 2259+586, in 
particular on shorter timescales than the exponential glitch 
relaxation timescale. This must be kept in mind when one 
considers the coincidence in timescales of the IE 2259+586 
rotational exponential decay and the average pulse profile’s 
relaxation to its preoutburst morphology. 

4.2.3. Comparison to SGRs 

From pulse frequency measurements leading up to and 
following the August 27 flare, it was determined that SGR 
1900+14 likely underwent rapid spin- down {Avfv~— lx 
10 -4 ) at or near the time of the flare (Woods et al. 1999b). An 
analysis of the pulsations during the flare itself (Palmer 2002) 
confirmed the existence of a timing anomaly, but only an 
upper limit on the timescale for the frequency change was 
measured (^1 day). For IE 2259+586, we detect a sudden 
increase in the rotation rate, as opposed to a decrease, in ad- 
dition to a much smaller magnitude than inferred for SGR 
1900+14. In general, the available limits on glitches (of either 
sign) in SGRs are poor. This is a result of fewer and less 
extensive phase-coherent timing solutions for these sources 
because of their strong timing noise (Woods et al. 2002) and 
weaker pulsed signals relative to most AXPs. Given the ob- 
servational limitations for monitoring glitch activity in the 
SGRs, it is not clear whether the glitch that accompanied the 
outburst in IE 2259+586 is common to SGR outbursts. 

Although there are not any lengthy phase-coherent solutions 
bridging the boundary between quiescent and burst-active 
states in SGRs to allow for identification of glitches, there are 
sufficient observations during and following a few outbursts to 
search for transient changes in spin-down rate. During the first 
2 months following the outburst of IE 2259+586, the frequency 
derivative may have changed sign (assuming the pulse shape 
changes did not affect our phase measurements; see § 3.1) and 
then rapidly switched back to spin-down with a higher mag- 
nitude compared with preoutburst. The average spin-down rate 
during the first 2 months after the glitch in IE 2259+586 was 
approximately double the preoutburst value. Similar behavior 
has not been seen during or following SGR outbursts. The fairly 
well sampled frequency histories of two SGRs (1900+14 and 
1 806—20) show that there is no direct correlation between burst 
activity and enhancements in the spin-down rate during or 
shortly following the outburst (Woods et al. 2002). In fact, the 
spin-down rate of SGR 1900+14 in the ~40 days following the 
1998 August 27 flare was at its lowest (in magnitude) historical 
value (Woods et al. 1999b, 2002). 

4.2.4. Postglitch Relaxation 

In radio pulsars, postglitch relaxation is almost always ob- 
served and can be classified into two regimes: one in which 
the initial frequency jump relaxes exponentially relatively 
quickly (i.e., timescales of hours to weeks) and one in which 
the frequency jump does not relax and which is generally 
accompanied by a permanent or very slowly (timescales of 
years) relaxing increase in the magnitude of the frequency 
derivative. Many radio pulsar glitches show both such behaviors 
(e.g., Shemar & Lyne 1996; Wang et al. 2000). These are 
today generally interpreted in terms of vortex creep theory, in 
which the glitches are due to a sudden impinning of angular 
momentum vortices in the neutron star crustal superfluid 
(Anderson & Itoh 1975). The latter, in the model, rotates faster 
than the crust but in equilibrium is loosely coupled to it via 
pinning of vortices to crustal lattice sites. In this picture, the 
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vortices creep outward slowly because of thermal activation 
(Alpar et al. 1984a, 1984c), slowly transferring angular mo- 
mentum to the crust. In equilibrium, the crust, slowed down 
by the external torque of magnetic braking, shares an identical 
spin-down rate with the superfluid in spite of a nonzero an- 
gular velocity differential between them. Vortex creep theory 
identifies the two observational postglitch relaxation regimes 
with different pinning regions (Alpar, Cheng, & Pines 1989, 
f 1993): the exponential relaxation is associated with superfluid 
regions with temperatures high compared with pinning ener- 
gies, so that equilibrium can be established with only a small 
angular velocity lag. This linear regime relaxation is mathe- 
matically equivalent to the original two-component model for 
the crust-superfliiid interaction suggested by Baym, Pethick, 
& Pines (1969), which did not incorporate the more recent 
vortex pinning theory. When die temperature is low compared 
with the pinning energy, a large lag is needed to establish 
equilibrium; this is the nonlinear regime and may be respon- 
sible for the longer timescale relaxations of spin-down rate 
enhancements. 

For IE 2259+586, qualitatively, we detect both such classical 
glitch behaviors as well. However, quantitatively, the postglitch 
relaxation of IE 2259+586 is very different from that seen in 
radio pulsars. Specifically, the exponential decay time is longer 
than what is observed in most pulsar glitches, and the Q value, 
the fraction of the glitch that heals, is also large by pulsar 
standards, although neither value is extreme. What is extreme is 
the combination: the net effect is that at 1 month postglitch die 
pulsar spins down for a couple weeks with more than double its 
long-term preoutburst pin-down rate. In radio pulsars, the 
spin-down enhancement is typically a few percent, nearly 2 
orders of magnitude smaller. This is because the pin-down 
enhancement is due to a temporarily reduced moment of inertia, 
as a linear component of the crustal neutron superfluid re- 
establishes equilibrium, and this superfiuid has at most a few 
percent of the moment of inertia of the crust. As discussed by 
Kaspi et al. (2003), a spin-down enhancement of order unity 
would imply that half the stellar moment of inertia was effec- 
tively decoupled from the crust following the glitch, much more 
than ought to be present in the crustal superfiuid. Indeed the 
crust itself has a relatively small moment of inertia compared 
with the core, which is thought to be coupled to the crust on 
short timescales (Alpar, Langer, & Sauls 1984b; Alpar & Sauls 
1988). Hence, the postglitch relaxation of IE 2259+586 could 
imply a core decoupling. This would provide the best evidence 
yet for the existence of core superfluicl 

However, the observed postglitch relaxation can be under- 
stood in terms of conventional crustal superfiuid if it is 
assumed that the superfiuid or crust angular velocity lag 
temporarily reversed at the time of die glitch (as was noted by 
Alpar, Pines, & Cheng 1990 in the case of the 1988 Christmas 
glitch of the Vela pulsar). The standard two-component theory 
has the following equations of motion for the crust (including 
r everything coupled to it on short timescales) and superfiuid: 

J ^ 4f(ftsf “ ^c) xr 

I c \l c = Acx, (3) 

* T 


4Af = — 


4f(ft S f flc) 
r 


( 4 ) 


where is the external torque on the crust and r is the crast- 
superfluid energy dissipation rate. The dissipation rate can be 


estimated from the postglitch exponential recovery timescale. 
From these equations, it follows that the change in the crustal 
angular frequency derivative is given by 


A tl c = 


4f A(fl s f Q c ) 
h ? 


( 5 ) 


or 


A(fi sf - n c ) = (6) 

hi 

For r ~ 15 days, 7 sf /7 c ~ 0.01, and ASl c ~2nx 10~ 14 rad 
s“ 2 , we find A(t2sf — Sl c ) ~ 8 x 10~ 6 rad s -1 . Since A(Qsf — 
Q c ) = Afigf “ AO c and AQ C =3x 10~ 6 rad s -1 , we have 
that AQsf a 5 x 10“ 6 rad s" 1 . As the lag was initially small 
and O c ~ 0.9 rad s“ ! , this implies the superfiuid angular 
velocity changed by only 1 part in 200,000. However, this 
change was about twice the observed change in the crust 
Thus, at the time of the glitch in this picture the superfiuid was 
temporarily spun down twice as much as the crust spun up. 
The observed enhanced spin-down postglitch is then due to 
the crust transferring angular momentum back to the super- 
fluid, as the latter attempts to arrange its vortices so as to 
reestablish equilibrium. 

The overall activity in glitches is a measure of the moment 
of inertia of the internal flywheel, which gives up angular 
momentum during glitch events (Link, Epstein, & Lattimer 
1999). The activity parameter is obtained by summing over all 
glitches measured over a long time interval I If the frequency 
lag between the superfiuid and the star is sourced only by the 
external torque, then one has 


W^T- 


( 7 ) 


In the case of IE 2259+586, we estimate Avjv = 4 x 10~ 6 
over an interval of T — 10 yr (the last outburst of IE 
2259+586 being estimated to have occurred 10 years earlier; 
Iwasawa et al. 1992; see also § 4.3). The resulting activity 
parameter is ~0.1, several times larger than what is measured 
in glitching radio pulsars (Link et al. 1999). Either 7 S f is 
proportionately larger or file spin rate of the superfiuid is being 
reduced by some mechanism other than vortex line unpinning 
and thermal creep. For example, smooth deformations of the 
stellar crust, driven by magnetic stresses, can have the latter 
effect (Thompson et al. 2000). 

A variety of models have been proposed to explain the or- 
igin of glitches in radio pulsars. In most models, the rotational 
lag between the superfiuid and the crust is the principal source 
of free energy. The frequency lag may as a result decrease in 
magnitude during the glitch, but it is more difficult to reverse 
its sign. For example, Ruderman (1991) suggested that file 
tension of the superfiuid vortices would fracture the crust, 
thereby allowing an outward shift in the positions of the 
vortices. Link & Epstein (1996) noted that any deposition of 
heat (on the order of ^lO 41 — 10 42 ergs in the deep crust) would 
accelerate the creep rate of the pinned vortex lines, thereby 
causing a transient spin-down of die superfluid and a spin-up 
of the rest of the star. Our observation of enhanced X-ray 
emission from IE 2259+586 is, indeed, consistent with this 
level of heat deposition in the crust. But it appears that this 
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cannot be the entire mechanism responsible for the observed 
glitch, because the vortex creep rate goes to zero as the spin 
rate of the superfluid approaches that of the crustal lattice. 

In a slowly rotating neutron star, magnetic stresses can act 
on timescales much shorter than the external torque. Glitch- 
like events result either from sudden fractures of the crust 
(Thompson & Duncan 1993) or from more gradual plastic 
deformations during which the vortices remain pinned 
(Thompson et al. 2000). The lowest energy deformations of 
this type are torsional. Consider, for example, a twisting 
motion of a circular patch of the crust that is offset in azimuth 
from the rotation axis. Starting from a uniform distribution, 
more vortices are advected outward than inward (because of 
the curvature of the stellar surface). A twist on the order of 
10 -2 rad is enough to provide the required spin-down of the 
crustal superfluid in IE 2259+586 (eq. [12] of Thompson et al. 
2000). By forcing the superfluid to rotate more slowly than it 
otherwise would, this mechanism also provides an explanation 
for the unusually large glitch activity that we inferred for IE 
2259+586 (see also Heyl & Hemquist 1999). 

Finally, as noted by Kaspi et al. (2003), the spin-down rate 
enhancement could also be due to a magnetospheric restruc- 
turing that causes the external torque to approximately double. 
However, given that the spin-up itself is inescapably enabled 
by an internal transfer of angular momentum, it would have to 
be a surprising coincidence that the external torque change 
offsets, to within a factor of 4-5, the internally generated 
frequency jump. This of course does not preclude some 
magnetospheric restructuring; indeed the residual pulse profile 
change following the initial rapid flux decay and the post- 
outburst infrared enhancement (Kaspi et al. 2003) suggest that 
likely occurred. However, it seems unlikely that a magneto- 
spheric restructuring is the cause of the enhanced spin-down 
just postoutburst. 

4 . 2 . 5 . Postglitch Long-Term Rotation 

Our glitch fitting clearly measured a long-term postglitch 
change in the spin-down rate. Interestingly, unlike that seen in 
a number of radio pulsars to date, the magnitude of the spin- 
down torque decreased. (In the case of the Crab pulsar, for 
example, file cumulative effect of the torque changes observed 
following some glitches is to increase the rate of spin-down by 
0.07% over 23 yr of observations; Lyne et al. 1993). If this is 
due to a superfluid effect, the conventionally observed long- 
term increase in magnitude of Cl c represents a reduction in the 
moment of inertia of the part of the neutron star that is coupled 
to its crust. This suggests that some fraction of the crustal 
superfluid is tightly pinned and not creeping outward and 
hence is decoupled from the stellar spin-down (Alpar et al. 
1994), with this fraction increasing at the glitch. The opposite 
must therefore be true of IE 2259+586 in this scenario. It has 
alternatively been suggested that the effect is due to a secular 
change in the magnetic moment (Ruderman 1991) or to a 
change in the electromagnetic torque acting on the star, driven 
by a reorientation of the magnetic moment with respect to the 
angular velocity of the star (Link & Epstein 1997). 

We now comment on the relative merits of a change in 
internal versus external torque in the case of IE 2259+586. 
The sign of the torque change is interesting, because there is 
independent evidence that the spin-down torque of IE 
2259+586 has decayed over time: its characteristic age, 
t c = 2 x 10 5 yr, is an order of magnitude greater than that of 
the SNR CTB 109 (estimated radius 12 pc; Kothes et al. 


2002) . A change in vortex creep rate with a timescale r much 
longer than 1 yr but much less than the characteristic age 
r c — 2 x 10 5 yr would appear as a permanent change in 
braking torque. The required moment of inertia of the super- 
fluid component would be 7 s f// = 0.05(r/r c ), comfortably 
below the fraction inferred for glitching radio pulsars (Link 
et al. 1999) unless the creep timescale were more than a 
million times longer than the observed 15 day relaxation. The 
simplest explanation for the required change is that the ex- 
ternal torque — a ~2.5% reduction in the magnetic moment — 
can be discounted. If, for example, the external magnetic field 
of IE 2259+586 were a centered dipole, the total field energy 
outside the star would amount to some 10 45 ergs, and the 
reduction in field energy would be a few hundred times larger 
that the excess X-ray emission we have observed. 

Nonetheless, changes in external torque can result from 
more subtle effects. If the magnetosphere is globally twisted 
because of the action of internal magnetic stresses, the 
resulting expansion of the poloidal field lines increases the 
current crossing the speed-of-light cylinder and therefore 
yields a greater torque than would exist from a simple dipolar 
field rotating in a vacuum (Thompson et al. 2002). Much 
of this effect results from the closed field lines that extend 
far from the star to a large radius R max , are anchored in a 
small fraction ~^(J?max/^Ns) -1 of its surface area, and contain 
only a small fraction ~(Rmax/^Ns) -3 of the energy in the 
nonpotential (toroidal) field outside the star. A slight relaxa- 
tion in the twist of the closed field lines is therefore consistent 
with the observed energetic output following the glitch if that 
reconfiguration occurs close to the magnetic poles. Persistent 
seismic activity in a magnetar will also modulate the spin- 
down torque, through a coupling of crustal shear waves to 
the magnetosphere, which in turns drives a particle wind 
(Thompson & Blaes 1998; Harding, Contopoulos, & Kazanas 
1999; Thompson et al. 2000). In this case, the observed 
change in torque corresponds to a very small reduction in the 
particle luminosity, by only ~ 10 31 ergs s” 1 . 

Finally, a change in torque could result from the excitation of 
precession during the glitch, followed by a permanent reori- 
entation of the magnetic axis relative to the rotation axis of the 
star (Link & Epstein 1997). That would require a precession 
angle of at least a few degrees, which is uncomfortably large if 
the nonspherical shape of the star is predominantly due to its 
magnetic field. The fractional difference in the sizes of the 
principal moments is ~10~ 6 [(5 2 )/(10 15 G) 2 ] (Wasserman 

2003) . However, the observed energy of the X-ray transient 
corresponds to less than ~10“ 3 of the total magnetic energy of 
the star (§ 4.2.1) and, thus, to a proportionately small change in 
the shape of the star. 

4.3. Implications for the Number of Active Magnetars 

As suggested by Kaspi et al. (2003), past reports of flux 
variability (Iwasawa et al. 1992; Baykal & Swank 1996), 
pulse profile changes (Iwasawa et al. 1992), and glitch activity 
(Heyl & Hemquist 1999) in IE 2259+586 likely indicate 
previous episodes of burst activity in this source. The best 
example of an inferred outburst from IE 2259+586 comes 
from a pair of Ginga observations in 1989 and 1990. Between 
1989 and 1990, the pulse profile of IE 2259+586 changed 
drastically, the flux increased by more than a factor 2, and the 
measured frequency showed evidence of a glitch. The close 
similarities with the well-studied outburst of 2002 June 
strongly suggests that a bursting episode from IE 2259+586 
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preceded the 1990 Ginga observation by roughly 1 week. 
Similar reports of flux variability and timing anomalies in the 
AXP be 1048.1—5937 (e.g., Oosteibroek et al. 1998) suggest 
past episodes of burst activity in this source as well. 

What is most intriguing about the detected and inferred 
outbursts in AXPs is that none have been detected with large 
FOV gamma-ray detectors (e.g., BATSE, Konus, Ulysses , 
etc.) that traditionally detect burst-active episodes in SGRs. 
This shows that we are missing low-intensity SGR-like out- 
bursts from magnetar candidates in our Galaxy (i.e., the SGRs 
are a sensitivity-limited sample). 

The number of active SGRs in our Galaxy was calculated 
previously to be ~10 (Kouveliotou et al. 1993). A key as- 
sumption in this estimation was that SGR outbursts are easily 
detected (i.e., there are always bright bursts within a given 
outburst that would trigger at least some of the large FOV 
gamma-ray instruments). The discovery of low-intensity SGR- 
like outbursts in AXPs clearly invalidates that assumption, 
making the estimate of 10 a lower limit only 

The future prospects for detecting weak SGR-like outbursts 
and addressing the question of the number of active SGRs in 
our Galaxy is promising. Until 2000 June, BATSE was die 
most sensitive large FOV detector for SGR burst emissions. 
The impending Swift mission and its sensitive BAT detector 
will be ~20 times more sensitive than BATSE was (D. Band 
2003, private communication). Just as the RXTE PCA opened 
our eyes to the larger population of low-luminosity SGR 
bursts (e.g., Gogu§ et al. 1999), perhaps Swift will reveal a 
previously unknown source population of dimm er SGRs. 

5. SUMMARY 

In 2002 June, the AXP IE 2259+586 was observed to emit 
more than 80 SGR-like bursts. Accompanying this outburst 
were several changes in the persistent X-ray emission prop- 
erties of the source, many of which are similar to what has 
been seen in SGRs, thus further blurring the distinction be- 
tween the two classes. We have quantified the observed 
changes by using data obtained with XMM-Newton and RXTE. 
In particular, we found the following: 

1. The flux increased by more than an order of magnitude 
and showed two decay components. The first component 
decayed rapidly within the first day of the outburst, while the 
second decayed much more slowly as a power law in time 
according to t~ 022 . 

2. The X-ray spectrum hardened during the outburst but 
almost fully recovered within 3 days. The spectrum at 21 days 
past the burst activity was significantly harder than at 7 days 
preoutburst but was fully consistent with historical spectral 
measurements of this source in quiescence. 

3. The phase dependence of the energy spectrum changed 
from that before to that after the outburst One week prior to the 
outburst, we observe significant variability in the photon index 
but not the blackbody temperature. Three days following the 
outburst, the dependence of the photon index on pulse phase 
flattened significantly. 
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4. The pulse profile changed suddenly during the observa- 
tion containing the burst activity during which much of the 
power moved to the fundamental frequency. The pulse profile 
rapidly returned to near its preoutburst shape within 1 week 
and showed only vay slow recovery thereafter. As with the 
flux and spectral changes, the recovery was not complete even 
after 1 year. 

5. The pulsed fraction decreased during the outburst to 
~15% but quickly recovered to the preoutburst value of ~23% 
within 6 days. 

6. IE 2259+586 suffered a glitch having an ordinary am- 
plitude [Ai/max/y = (4.24 ± 0.11) x 10“*] but a unique re- 
covery. Approximately 19% of the glitch recovered on a 
timescale of ~16 days, although the recovery was not exactly 
exponential in form. During the recovery, the torque was en- 
hanced relative to the preoutburst value by a factor of ~2. We 
detect a long-term ~2% reduction in the spin-down rate fol- 
lowing the glitch [A z> = (2.18 ± 0.25) x 10~ 16 Hz s -1 ]. 

7. The measured glitch epoch precedes the observed burst 
activity by ~12 hr. Given the rapid flux decay during the 
outburst, die true onset of this burst activity may have followed 
the glitch. 

The cumulative properties of the outburst in IE 2259+586 
lead us to conclude that the star suffered some major event that 
was extended in time and had two components, one tightly 
localized on the surface of the star (i.e., a fracture or a series of 
fractures) and the second more broadly distributed (possibly 
involving a smoother plastic change). This event affected both 
the superfluid interior and the magnetosphere. The glitch 
points toward a disturbance within the superfluid interior, 
while the extended flux enhancement and pulse profile change 
suggest an excitation of magnetospheric currents and crustal 
heating. Finally, we show that the lack of detection of AXP 
outbursts with all-sky gamma-ray detectors suggests there 
exists a larger population of active SGRs in our Galaxy than 
previously thought 
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